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Abstract 

This work is devoted to the study of the initial boundary value problem for a general non 
isothermal model of capillary fluids derived by J.E Dunn and J.Serrin (1985) in |X5|, I25j. 
which can be used as a phase transition model. 

We distinguish two cases, when the physical coefficients depend only on the density, and 
the general case. In the first case we can work in critical scaling spaces, and we prove 
global existence of solution and uniqueness for data close to a stable equilibrium. For 
general data, existence and uniqueness is stated on a short time interval. 
In the general case with physical coefficients depending on density and on temperature, 
additional regularity is required to control the temperature in L°° norm. We prove global 
existence of solution close to a stable equilibrium and local in time existence of solution 
with more general data. Uniqueness is also obtained. 

1 Introduction 

1.1 Derivation of the Korteweg model 

We are concerned with compressible fluids endowed with internal capillarity. The model 
we consider originates from the XlXth century work by van der Waals and Korteweg 
|21j and was actually derived in its modern form in the 1980s using the second gradient 
theory, see for instance [20l [26] . 

Korteweg-type models are based on an extended version of nonequilibrium thermodynam- 
ics, which assumes that the energy of the fluid not only depends on standard variables 
but on the gradient of the density. Let us consider a fluid of density p > 0, velocity 
field u G M N (N > 2), entropy density e, and temperature 9 = (|§ ) p . We note w = Vp, 
and we suppose that the intern specific energy, e depends on the density p, on the en- 
tropy specific s, and on w. In terms of the free energy, this principle takes the form of a 
generalized Gibbs relation : 

de = Tds + -^dp H — 6* ■ dw 
p l p 
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where T is the temperature, p the pressure, cf> a vector column of R and 4>* the adjoint 
vector. 

In the same way we can write a differential equation for the intern energy per unit volume, 
E = pe, 

dE = TdS + gdp + <fi* ■ dw 

where S = ps is the entropy per unit volume and g = e — sT+ ^ is the chemical potential. 
In terms of the free energy, the Gibbs principle gives us: 

dF = -Sdf + gdp + (/)* ■ dw. 

In the present chapter, we shall make the hypothesis that: 

(j) = KW. 

The nonnegative coefficient k is called the capillarity and may depend on both p and T. 
All the thermodynamic quantities are sum of their classic version (it means independent 
of w) and of one term in \w\ 2 . 

In this case the free energy F decomposes into a standard part Fq and an additional term 
due to gradients of density: 

F = F + -k\w\ 2 . 

We denote v = - the specific volume and k = vk. Similar decompositions hold for S, p 
and g: 

P = Po~ 2 K p\ w \ 2 wnere: k p = k 'v and Po = -{fo)' v 

9 = 9o + 2 K 9\ W \ 2 wnere: K e = k- fk' f and e = f - T(/ )~. 

The model deriving from a Cahn-Hilliard like free energy (see the pioneering work by 
J.E.Dunn and J.Serrin in [15j and also in [HEKTT]), the conservation of mass, momentum 
and energy read: 

' d t p + div(pu) = 0, 

< dt(pu) + div(/ni (g> u +pl) = div(K + D) + pf, 

k d t (p(e + \u 2 )) + divKpe + \p\u\ 2 + p)) = div((£> + K) ■ u - Q + W) + pf ■ u, 
with: 

D = (Adivn)/ + p(du + Vu), is the diffusion tensor 
K = {pd\v(j))I — (fiw*, is the Korteweg tensor 
Q = — ryVT, is the heat flux. 

The term 

W = (d t p + u* ■ Vp)(j) = -(pdivu)<j> 

is the intersticial work which is needed in order to ensure the entropy balance and was 
first introduced by Dunn and Serrin in |15j . 

The coefficients (A, p) represent the viscosity of the fluid and may depend on both the 
density p and the temperature T. The thermal coefficient r\ is a given non negative 
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function of the temperature T and of the density p. 

Differentiating formally the equation of conservation of the mass, we obtain a law of 
conservation for w. 

dfW + div(uw* + pdu) = . 

One may obtain an equation for e by using the mass and momentum conservation laws 
and the relations: 

div((-p7 + K + D)u) = (div(-p7 + K + D))-u- pd\w{u) + (K + D) : Vu . 
Multiplying the momentum equation by u yields: 

(div(-pJ + K + D))-u= (dt(pu) + divfrwm*)) • u = d t {^-) + dw(^u) . 
We obtain then: 

p(d t e + u* ■ Ve) +pdivn = (K + D) : Vu + div(W - Q) . 
In substituting K, we have (with the summation convention over repeated indices): 

K : Vu = pdiv^divu — 4>iWjdiUj , 

while: 

—divW = div((pdivn)0) = /)(div(/>)(dhat) + (w* ■ <^>)divu + pdj^Uj. 

In using Wj = djp, we obtain: 

K : Vu — dwW = —(w* ■ </>)divu — (pidj(pdiUj) 

= — (w* • (j))divu — (div(pdu)) • <fi. 

Finally, the equation for e rewrites: 

p(dte + u* • Ve) + {jp + w* • </>)divn = D : Vu — (div(pdu)) • <fi — divQ . 
From now on, we shall denote: dt = dt + u* ■ V. 

1.2 The case of a generalized Van der Waals law 

From now on, we assume that there exist two functions Ho and 111 such that: 

P0 = fu[(v)+U' (v), 

eo = -IIo(v) + <p(T)-T<p'(T). 

We now suppose that the coefficients A, p depend on the density and on the temperature, 
and in all the sequel the capillarity k doesn't depend on the temperature. 
Moreover we suppose that the intern specific energy is an increasing function of T: 

(A) *'(T) >0 with *(T) = 99(f) -Ttp'(T). 

We then set 6 = *(T) and we search to obtain an equation on 9. In what follows, we 
assume that k depends only on the specific volume. 



Obtaining an equation for 9 : 

As: 



we thus have: 



e = -U (v) + 8 + ^k\w\ 2 , 
dte = —U (v)d t v + d±0 + -K v \w\ 2 dtv + kw* ■ dtw . 



By a direct calculus we find: 

dtv = vdivu and w* ■ dtw = — \w\ 2 divu — div(pdu) ■ w . 

Then we have: 

dfO = dte + v(p — TII 1 (t;))divti + K|u>| 2 divu + ndiv(pdu) • w. 
And in using the third equation of the system, we get an equation on 9: 

dtO + vdivQ + wTIT 1 (t;)divu = vD : Vu + div(pdu) • (kw — v(f>) + (k\w\ 2 — vw* • <p)divu . 
But as we have 4> = K w and k = v k we conclude that: 

d t 9 - vdw( X V9) + v^' 1 (9)u' 1 (v)dWu = vD : Vu 
with: x(p,0) = v(P,T)(^- 1 )'(9) . 

Obtaining a system for p, u, 9: 

We obtain then for the momentum equation: 

divD V Po divK lV(K p \w\ 2 ) 
d t u 1 = h ■ 



P P p 2 p 

where K p = k — pn' p . 

And by a calculus we check that: 

divK + \v(K p \w\ 2 ) = P^(kA p ) + ^V(K p \Vp\ 2 ). 

Indeed we have: 

I =V(pdiv(/cVp)) - div(Kww*) + ^V(K p \w\ 2 ) 

= [pV(KAp) + KVpAp + pV(K p \p\ 2 ) + K p \Vp\ 2 Vp] - [k<Hv(ww*) + K p w Vpw] 

I 

+ [^VM 2 + yVl^l'Vp - ^V( P K p \w\ 2 )], 



= [pV(«Ap) + KVpAp+ ^V(k' p \p\ 2 )} - [/ediv(W)] + [|v|w| 2 ], 



--[pV(KAp) + KVpAp + ^V(^|p| 2 )] - [Kwdivw + |v|u>| 2 ] + [^V|u>| 2 ], 

-pV(kAp) + |v( k ;|vh 2 ) • 
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Finally we have obtained the following system: 
d t p + div(pu) = 0, 



(NHV) 



= V(^|V P | 2 ), 

dt e + u .ve- div (* w > + W ^MdivH = 

p p p 



where: P = U' Q , Pi = n' x and T = ^-\6). 
We supplement (NHV) with initial conditions: 

p t =o = Po > u t=0 = n , and t=o = O . 
1.3 Classical a priori-estimates 

Before getting into the heart of mathematical results, let us derive the physical energy 
bounds of the (NHV) system when n is a constant and where the pressure just depends 
on the density to simplify. Let p > be a constant reference density, and let II be defined 
by: 



U(s) 



Jo Z 



so that P (s) = sU (s) - U(s) , U (p) = and: 

d t YL(p) + div(«n(p)) + P (p)div(n) = in V((0, T) x R N ). 

Notice that II is convex as far as P is non decreasing (since Pq(s) = sll" (s)), which is the 
case for 7-type pressure laws or for Van der Waals law above the critical temperature. 
Multiplying the equation of momentum conservation in the system (NHV) by pu and 
integrating by parts over R^, we obtain the following estimate: 

/ (\p\u\ 2 + pe + (Ii(p)-Ii(p)) + ^\Vp\ 2 )(t)dx + 2 f I (2fiD(u) : D(u) 
Jr n z z jo jr n 

+ (A + p)\divu\ 2 )dx < [ + P0 9 + (U( P0 ) - U(p)) + ^|Vpo| 2 )^. 

J r n Zp Z 

It follows that assuming that the initial total energy is finite: 

|2 



eo 



m 



9 + Po e + (u(p ) - n(p)) + -\v Po \ 2 )dx < +cx) , 



then we have the a priori bounds: 

U(p)-U(p), p\u\ 2 , and p6 G L°°(0, 00, L 1 (M^)) 
VpG L co (0,oo,L 2 (R N )) N , and Vu G L 2 (0, 00, R N ) N \ 
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2 Mathematical results 

We wish to prove existence and uniqueness results for {NHV) in functions spaces very 
close to energy spaces. In the non isothermal non capillary case and P(p) = ap" 1 , with 
a > and 7 > 1, P-L. Lions in [22] proved the global existence of variational solutions 
(p, u, 6) to (NHV) with k = for 7 > f if N > 4, 7 > ^ if N = 2 > 3 and initial data 
(po,mo) such that: 



These solutions are weak solutions in the classical sense for the equation of mass conser- 
vation and for the equation of the momentum. 

On the other hand, the weak solution satisfies only an inequality for the thermal energy 
equation. 

Notice that the main difficulty for proving Lions' theorem consists in exhibiting strong 
compactness properties of the density p in L p loc spaces required to pass to the limit in the 
pressure term P(p) = ap 1 . 

Let us mention that Feireisl in [16] generalized the result to 7 > y in establishing that we 
can obtain renormalized solution without imposing that p £ Lf , for this he introduces 
the concept of oscillation defect measure evaluating the lost of compactness. 
We can finally cite a very interesting result from Bresch-Desjardins in [5], [6] where they 
show the existence of global weak solution for (NHV) with k = in choosing specific 
type of viscosity where p and A are linked. It allows them to get good estimate on the 
density in using energy inequality and to can treat by compactness all the delicate terms. 
This result is very new because the energy equation is verified really in distribution sense. 
In [23], Mellet and Vasseur improve the results of Bresch,Desjardins in generalize to some 
coefficient p and A admitting the vacuum in the case of Navier-Stokes isothermal, they 
use essentially a gain of integrability on the velocity. 

In the case k > 0, we remark then that the density belongs to L°°(0, 00, ii" 1 (IR^)). Hence, 
in contrast to the non capillary case one can easily pass to the limit in the pressure term. 
However let us emphasize at this point that the above a priori bounds do not provide any 
L°° control on the density from below or from above. Indeed, even in dimension N = 2, 
H 1 functions are not necessarily locally bounded. Thus, vacuum patches are likely to 
form in the fluid in spite of the presence of capillary forces, which are expected to smooth 
out the density. Dan chin and Desjardins show in [H] that the isothermal model has weak 
solutions if there exists c\ and M\ such that: 



The vacuum is one of the main difficulties to get weak solutions, and the problem remains 
open. 

In the isothermal capillary case with specific type of viscosity and capillarity p(p) = pp 
and X(p) = 0, Bresch, Desjardins and Lin in [7] obtain the global existence of weak 
solutions without smallness assumption on the data. We can precise the space of test 
functions used depends on the solution itself which are on the form pcj) with G C^(R N ). 




ci < \p\ < Mi and \p - 1| << 1. 



6 



The specificity of the viscosity allows to get a gain of one derivative on the density: 
pi /.-(// 2 ). ' 



Existence of strong solution with k, p and A constant is known since the work by Hat- 
tori an Li in [18], [19] in the whole space K^. In [14] . Danchin and Desjardins study 
the well-posedness of the problem for the isothermal case with constant coefficients in 
critical Besov spaces. 

Here we want to investigate the well-posedness of the full non isothermal problem in 
critical spaces, that is, in spaces which are invariant by the scaling of Korteweg's system. 
Recall that such an approach is now classical for incompressible Navier-Stokes equation 
and yields local well-posedness (or global well-posedness for small data) in spaces with 
minimal regularity. 

Let us explain precisely the scaling of Korteweg's system. We can easily check that, if 
(p,u,9) solves (NHV), so does (p\,u\,9\), where: 

p\(t,x) = p(X 2 t, Xx) , u\(t,x) = Xu(X 2 t, Xx) and 9\(t,x) = X 2 9(X 2 t, Xx) 

provided the pressure laws Poj Pi have been changed into A 2 i-b> A 2 Pi. 

Definition 2.1 We say that a functional space is critical with respect to the scaling of 
the equation if the associated norm is invariant under the transformation: 

(p,u,6) — ► (p\,u x ,9 x ) 

(up to a constant independent of X). 

This suggests us to choose initial data (po,uo,9o) i n spaces whose norm is invariant for 
all A > by (po,t*o,0o) — (Po(A-), Au (A-), X 2 9 (X-)). 

A natural candidate is the homogeneous Sobolev space 

but since H N / 2 is not included in L°°, we cannot expect to get L°° control on the density 
when po G H N/2 . The same problem occurs in the equation for the temperature when 
dealing with the non linear term involving , I'~ 1 (#). 

This is the reason why, instead of the classical homogeneous Sobolev space we 
will consider homogeneous Besov spaces with the same derivative index B s = .BliO^O 
(for the corresponding definition we refer to section 4). 

One of the nice property of B s spaces for critical exponent s is that B N / 2 is an algebra 
embedded in L°°. This allows to control the density from below and from above, without 

JV 

requiring more regularity on derivatives of p. For similar reasons, we shall take 9q in B 2 
in the general case where appear non-linear terms in function of the temperature. 
Since a global in time approach does not seem to be accessible for general data, we will 
mainly consider the global well-posedness problem for initial data close enough to stable 
equilibria (Section 5). This motivates the following definition: 

Definition 2.2 Let p > 0, 9 > 0. We will note in the sequel: 



One can now state the main results of the paper. 

The first three theorems concern the global existence and uniqueness of solution to the 
Korteweg's system with small initial data. In particulary the first two results concern 
Korteweg's system with coefficients depending only on the density and where the intern 
specific energy is a linear function of the temperature. 

Theorem 2.1 Let N > 3. Assume that the function ^ defined in (A) satisfies ^f(T) = 
AT with A > and that all the physical coefficients are smooth functions depending only 
on the density. Let p > be such that: 

K(p) > 0, p(p) > 0, X(p) + 2p(p) > 0, r](p) > and d p P (p) > 0. 

Moreover suppose that: 

~JV_i A N i _, ~N i ]V_n 

qo £ B z ^.MoeBi \ T eBi x, 2 

There exists an £q depending only on the physical coefficients (that we will precise later) 
such that if: 

then (NHV) has a unique global solution (p,u,T) in E N / 2 where E s is defined by: 

E s =[C b (R + ,B s - 1 ' s )nL 1 {R + ,B s+1 ' s+2 )} x [C b (R+, B 3 ' 1 ^ n L\R + , B S+1 ) N ] 
x [C b {R+, b 3 - 1 ' 3 - 2 ) n L\R+, B 3+1 ' 3 )}. 

Remark 1 Above, B 3,t stands for a Besov space with regularity B 3 in low frequencies 
and B 1 in high frequencies (see definition \3. 3\) . 

The case N = 2 requires more regular initial data because of technical problems involving 
some nonlinear terms in the temperature equation. 

Theorem 2.2 Let N = 2. Under the assumption of the theorem \2.1\ for ^ and the 

physical coefficients, let e > and suppose that: 

q E B°' 1+£ ' , u & B°' £ ' , T E B°'- 1+£ ' . 
There exists an eo depending only on the physical coefficients such that if: 

Ikollgo.i+e' + ll n o||^o, £ ' + ll^ollgo.-iW - £ o 
then (NHV) has a unique global solution (p,u,T) in the space: 

E' =[C b (R + ,B ' 1+£ ')r\L 1 (R + ,B 2 ' 3+£ ')} x [C b {R + ,B ' £ ') 2 r\L 1 (R + ,B 2 ' 2+£ ') 2 ] 
x [C b {R+ , B^- 1+£ ) n L 1 (M+ , B 2 ' 1+£ ' )}. 

In the following theorem we are interested by showing the global existence of solution for 
Korteweg's system with general conditions and small initial data. In order to control the 
non linear terms in temperature more regularity is required. That's why we want control 
the temperature in norm L°°. 
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Theorem 2.3 Let N > 2. Assume that ^ be a regular function depending on 9. Assume 
that all the coefficients are smooth functions of p and 9 except k which depends only on 
the density. Take (p,T) such that: 

K {p) > 0, p{p, f) > 0, X(p, f ) + 2p{p, f) > 0, n{p, f) > and d p P {p, f) > 0. 

Moreover suppose that: 

~N 1 JV ,1 ~JV -I JV ~JV i JV 

q G B~z + , u G B-' 1 '-, T G BT -1 'T. 
There exists an e\ depending only on the physical coefficients such that if: 

On ~jv i jv i-r + ktn ~jv -, jv + Xn ~jv , jv < £ 
f/ien (NHV) has a unique global solution (p,u,T) in: 

JV . „ . ~JViJV, 1n i , ~ JV i 1 JV i q. , . „ . ~K i W, if 

F~ =[Ci,(lR + , B 2 n L 1 (K + , B^ +1 -- +3 )] x [C&(]R+, B'a" '"a") 
nL 1 (M + ,Bf +1 'f +2 )^][a(R + , J Bf- 1 'f)nL 1 (IR + , J Bf +1 'f +2 )]. 

In the previous theorem we can observe for the case N = 2 that the initial data are very 
close from the energy space of Bresch, Desjardins and Lin in [7]. 

In the following three theorems we are interested by the existence and uniqueness of 
solution in finite time for large data. We distinguish always the differents cases N > 3 
and N = 2 if the coefficients depend only on p, and the case where the coefficients depend 
also on T. 

Theorem 2.4 Let N > 3, and ^ and the physical coefficients be as in theorem \2.l\ We 
suppose that (qo,uo,To) G B~ x (B~~ l ) N x B~~ 2 and that po > c for some c > 0. 
Then there exists a time T such that system (NHV) has a unique solution in Ft 

F t ={C t (B%) n 4(flf + 2 )] x [CriB?- 1 ) 11 n L T {B% +1 ) N ] 
x [C r (Bf" 2 )n4(BT)] . 

For the same reasons as previously in the case N = 2 we can not reach the critical level 
of regularity. 

Theorem 2.5 Let N = 2 and e > 0. Under the assumptions of theorem \2.1\ for ^ and 

the physical coefficients we suppose that (qo,uo,To) G B 1,1+£ x (B°' £ ) 2 x B >~ 1+£ and 
Po > c for some c > 0. 

Then there exists a time T such that the system has a unique solution in -Ft (2) with: 

F T (2) =[C T {B^ l+£, )r\L T {B^ +£ ')} x {C T (B > £ ') 2 r\L T (B 2 > 2+£ ') 2 ] 
x [C T {B-^- 1+£ ')r\LT{B^ 1+£ ')} . 

In the last theorem we see the general system without conditions, and like previously we 
need more regular initial data. 
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Theorem 2.6 Under the hypotheses of theorem \2.3\ we suppose that: 

~N_ N_ at JV 

(qo, uq, T ) G B 2 > 2 x (B 2 ) JV x 2 and po > c for some c > 0. 
T/ien £/iere exists a time T such that the system has a unique solution in : 

F' T =[c T (m^ +l )^L l T {m +2 ^)} x [c t (bt fmUs^Y] 

x [C T (B% ) r\L l T {B^ +2 )] . 

This chapter is structured in the following way, first of all we recall in the section [3] 
some definitions on Besov spaces and some useful theorem concerning Besov spaces. 
Next we will concentrate in the section 0] on the global existence and uniqueness of 
solution for our system (NHV) with small initial data. In subsection 14.11 we will give 
some necessary conditions to get the stability of the linear part associated to the system 
(NHV). In subsection 14.21 we will study the case where the specific intern energy is linear 
and where the physical coefficients are independent of the temperature. In our proof we 
will distinguish the case N > 3 and the case N = 2 for some technical reasons. In the 
section [5] we will examine the local existence and uniqueness of solution with general 
initial data. For the same reasons as section H] we will distinguish the cases in function 
of the behavior of the coefficients and of the intern specific energy. 

3 Littlewood-Paley theory and Besov spaces 
3.1 Littlewood-Paley decomposition 

Littlewood-Paley decomposition corresponds to a dyadic decomposition of the space in 
Fourier variables. 

We can use for instance any ip G C°°(R N ), supported in C = {£ G R N /j < |£| < §} such 
that: 

5>(2-'o = i if e/o. 

Denoting h = T~ x <p, we then define the dyadic blocks by: 

Aiu = 2 lN / h(2 l y)u(x - y)dy and Siu = A k u . 

Formally, one can write that: 

u = ^2 A k u . 

This decomposition is called homogeneous Littlewood-Paley decomposition. Let us ob- 
serve that the above formal equality does not hold in S (M. N ) for two reasons: 

1. The right hand-side does not necessarily converge in S' (M N ). 

2. Even if it does, the equality is not always true in S' (M N ) (consider the case of the 
polynomials). 

However, this equality holds true modulo polynomials hence homogeneous Besov spaces 
will be defined modulo the polynomials, according to [1]. 
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3.2 Homogeneous Besov spaces and first properties 
Definition 3.3 For s£l, and u £ S'(R N ) we set: 

1Mb' = (^( 2 '1 A HIlO- 

lez 

A difficulty due to the choice of homogeneous spaces arises at this point. Indeed, ||.||b s 
cannot be a norm on {u £ S (R N ), \\u\\b s < +00} because ||u||b s = means that u is 
a polynomial. This enforces us to adopt the following definition for homogeneous Besov 
spaces, see [I]. 

Definition 3.4 Let s £ E. 

Denote m = [s — y] i/ s — y ^ ^ anc ^ m = s ~~ y ~~ 1 otherwise. 

• If m < 0, then we define B s as: 

B s = lue S'(R N ) /\\u\\ B > < 00 and u = ^ Am inS'(R N ) i . 

• If rn > 0, we denote by PmjR^] i/ie set of polynomials of degree less than or equal 
to m and we set: 

B s = lueS'(R N )/V m [R N ] /\\u\\ B » <oo and u = A t u in S' {R N )/V m [R N } \ . 

lei ' 

Proposition 3.1 The following properties hold: 

1. Density: If \s\ < y > then is dense in B s . 

2. Derivatives: There exists a universal constant C such that: 

C _1 |l u ll-B s ^ II V^IIb 8 - 1 < C|M|.B S - 

3. Algebraic properties: For s > 0, B s n L°° is an algebra. 

4. Interpolation: (B Sl ,B S2 ) d l = B ds ' 2+ ^ d ^ Sl . 

3.3 Hybrid Besov spaces and Chemin-Lerner spaces 

Hybrid Besov spaces are functional spaces where regularity assumptions are different in 
low frequency and high frequency, see |12j . 
They may be defined as follows: 

Definition 3.5 Let s, t £ R. We set: 

q<0 q>0 

Let m = — + 1 — s], we then define: 
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• B s > 1 = {u£ S'(R N ) /\\u\\ Sstt < +00}, ifm < 

• = {u e S'{R N )/V m [R N ] /\\u\\ Ss!t < +00} ifm > 0. 

Let us now give some properties of these hybrid spaces and some results on how they 
behave with respect to the product. The following results come directly from the parad- 
ifferential calculus. 

Proposition 3.2 We recall some inclusion: 

• We have B s > s = B s . 

• Ifs<t then B s ' f = B s f] B t , if s > t then B 3 ^ = B s + B f . 

• If S! < s 2 and h > t 2 then B S1 ^ ^ B S2 * 2 . 

Proposition 3.3 For all s, t > 0, we have: 

\\uv\\^ tt < C(\\u\\ L °°\\v\\g S:t + |MU°°Nlj?<>>t) . 
For all s\, s 2 , t\, t 2 < y such that min(si + s 2 , t\ + 1 2 ) > we have: 

\\uv\\ ~„ , t N „ , , JV < C\\u\\ 5 S1 u \\V II s s „ to . 

For a proof of this proposition see [T2] • We are now going to define the spaces of Chemin- 
Lerner in which we will work, which are a refinement of the spaces: 

L P T {B S ) := L p (0,T,B s ). 



Definition 3.6 Let p € [1, +00 [, T £ [1, +00] and s £ E. We then denote: 

hWlMSn*) =E 2 ^( f \\^n{t)\\ P L2 dt) l/p + ^ lS2 { f \\^lu{t)\\" L2 dt) l/p . 
l<0 Jo l>0 Jo 

And we have in the case p = 00: 

K0 l>0 

We note that thanks to Minkowsky inequality we have: 

IHIz4,(5 a i> s 2) ^ \\ u \\l p t (B s ±< s 2) and W U W L^B*!^) = W U Wl^(B s 1' s 2)- 
From now on, we will denote: 



\u\\~ 



LUB' 



1) = E 2tol (/ \\^(t)\\ p LP dty 

K0 J ° 
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<^ 2) = £ 2 ' S2 (/ \\^(t)\\^) 1/p - 

l>0 JU 



Hence: 



We then define the space: 

L p T (B s ^) = {u G Z4(S™)/||u|| £ , (Sn , S2) < 00} . 

We denote moreover by Ct{B Si,S2 ) the set of those functions of Lj?\B Sl,S2 ) which are 
continuous from [0, T] to B Sl,S2 . In the sequel we are going to give some properties of 
this spaces concerning the interpolation and their relationship with the heat equation. 

Proposition 3.4 Let s, t, s±, S2 G M., p, pi, P2 G [l,+oo]. We have: 

1. Interpolation: 

^II^B..')^^!!^^!.*!)!!"!!^^) «^0€[O,1] and - = _ + __, 

a = 0«i + (1 - 9)8 2 , t = 6t 1 + (l- 6)t 2 . 

2. Embedding: 



L P T {B S ^) ^ L p t (Cq) and C T (B^) ^ C([0,T] 



x 



The Lj,(Bp) spaces suit particulary well to the study of smoothing properties of the heat 
equation. In [9], J-Y. Chemin proved the following proposition: 



Proposition 3.5 Let p G [1, +00] and 1 < p% < p\ < +00. Let u be a solution of: 

dtu — pAu = f 
u t =o = u . 

Then there exists C > depending only on N, p, pi and p 2 such that: 

W U \\Z P1 (B s + 2 /Pl) — Cll^ollB 8 + C\\f\\j j P2^ B s-2+2/ P2 y 

To finish with, we explain how the product of functions behaves in the spaces of Chemin- 
Lerner. We have the following properties: 

Proposition 3.6 

Let s > 0, t> 0, l/p 2 + 1/P3 = 1/pi + 1/pa = 1/p < 1, u G L£ 3 (£? s '*) n L%{L°°) and 
v G Jj^(B s,t ) n (L°°). 

Then uv G Lj,(B s,t ) and we have: 

1 1 uv 1 1 l p T (§*>*) - C \\ u \\l p T l {l^)\\ v \\l p t a (b«>*) + H' u llz^(L°°)ll' u llL^(B s > t )' 
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Ifsi,s 2 MM < f , + > 0, h + t 2 > 0, i + i = i < 1, « G L P ^{B S ^) and 



v G L p2 (B s ^ t2 ) then uv G L^(B sl+S2 ~f' tl+t2 'f) and: 

For a proof of this proposition see [12] . Finally we need an estimate on the composition 
of functions in the spaces Lj,(Bt) (see the proof in the appendix). 

Proposition 3.7 Let s > , p G [l,+oo] and u l7 u 2 - • • , Ud G L^(Bp) n L^(L°°). 

(i) Let F G W^ 2 ' 00 ^) suc/t tfwrf F(0) = 0. T/ien F(«i, u 2 , ■••,««*)€ 

More precisely, there exists a constant C depending only on s,p,N and F such that: 

\\F(ui, U2, ■ ■ ■ ,Ud) ^ C(ll n l||L5?(L°°) 5 II -u 2||l5?(L°°) 5 • • • > \\ Ud \\ ) 



fiij Lei it G L P T (B S ^ S2 ), s 1} s 2 > i/ien we We F(u) G L^(B Sl ' S2 ) and 

(Hi) Ifv,u£ L p T (B s p ) n L^(L°°) and G G w/^ 3 ' 00 ^) i/ien G(u) - G(v) belongs to 
Lj,(Bp) and there exists a constant C depending only of s,p,N andG such that: 

\\G(u) - G(v)\\ Z P T{BSp) < CQMItoQ^oo), ||w||l«.(z,oq))(||u - «||^, (B .)(1 + 

+ \Hl ¥ (l°°)) + \\v - «IU|?(L«»)(l|tt|lz5,(fl.) + Hr§,(B-)))- 

If v, u € L p T (B s p uS2 ) n L~(L°°) and G G W^ 3 ' 00 ^) ffcen G(n) - 6e/on 5 s to 
L P -,(Bp 1 ' S2 ) and it exists a constant C depending only of s,p,N andG such that: 

\\G(u) - G(v)\\i^ S n,s^ < C([|u[| L oo(x,oo), ||u|| jL oo( L o 0) )(||w - ( _g»i.»2)(l + [MUf?(z,°°) 

+ IMU»(i<»)) + Ik - u \\l™(l™)(\\u\\zp t (b s /- 82 ) + IMIz^S^ 2 )))' 

The proof is an adaptation of a theorem by J.Y. Chemin and H. Bahouri in [2], see the 
proof in the Appendix. 

4 Existence of solutions for small initial data 
4.1 Study of the linear part 

This section is devoted to the study of the linearization of system (NHV) in order to get 
conditions for the existence of solution. We recall the system (NHV) in the case where 
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k depends only on the density p: 
' d t p + dW(pu) = 0, 
d t (pu) + divO <g> u) - divD - p V(kA p ) + (V(P (p) + TP^p)) 



(NHV) 



pV(^|vh 2 ), 



^_dMxW) +r ^) div(u) 



L> : Vn 



P P P 

Moreover we have: 

div(D) = (A + ^)Vdivu + pAu + V(A)divu + (du + Vu)Vp, 

= (2A + p)Vdivu + pAu + «9iA(p, 6»)Vpdivn + d 2 X(p, 6)V6divu 
+ (du + Vu)d 1 p(p, 9)X7p + (du + Vu)d 2 p(p, 0)V6. 

We transform the system to study it in the neighborhood of (p, 0, 9). Using the notation 
of definition 12.21 we obtain the following system where F, G, H contain the non linear 
part: 

' dtq + divu = F , 

d t u- ^Au- ( A + A) Vdiv u - pRVAq + (P (p) + TP (p))Vq 
P P 

+ ^VT = G , 
pip (T) 



(M) 



d t T-*AT + ^±dWu = H. 
P P 

This induces us to study the following linear system: 
dtq + divu = F 

d t u - pAu - (p, + A) - eVAq - (5Vq - 7VT = G 
d t T - aAT + Sdivu = H 
d t u - pAu = PG 



(M' 



where u, e, a, /3, 7, 5 and p are given real parameters. Note that system (M) with right 
hand side considered as source terms enters in the class of models (M ), it is only a 
matter of setting: 



P = ^, A = ^, e = pn, (3 = Pq(p) + TP[(p), 7 
P P 



Pi(p) = X <? 

P^(fy a p'° 



tpi(p) 



We transform the system in setting: 

d = A _1 divu and f2 = A~ 1 curlu 
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where we set: A s h = ^ 7 ~ 1 (|^| s /i) (the curl is defined in the appendix). 

We finally obtain the following system in projecting on divergence free vector fields and 

on potential vector fields: 

( d t q + Ad = F, 
d t d - v&d - eA 3 q - (5 Aq - 7 AT = A _1 div G, 



(M[) 



d t T - aAT + 5 Ad = H, 

d t n-jiAn = A _1 cuxiG, 

u = -A _1 Vd- A _1 divfi. 



The last equation is just a heat equation. Hence we are going to focus on the first three 
equations. However the last equation gives us an idea of which spaces we can work with. 
The first three equation can be read as follows: 



(M' 2 



mo \ / mo 

d t | d(t,o +A(0 d(t,o 



A _1 div G(t, 
H(t,0 



where we have: 



MO 



o iei o 

-e\e-m\ »\e — riei 

8\Z\ a\C\ 2 



The eigenvalues of the matrix — A(0 are of the form |^| 2 A^ with being the roots of 
the following polynomial: 

Pt(X) = X 3 + (u + a)X 2 + (e + va + 1 ^^) X + (« e + j^) ■ 

For very large £, the roots tend to those of the following polynomial (by virtue of conti- 
nuity of the roots in function of the coefficients): 

X 3 + (v + a)X 2 + (e + ua)X + ae. 



The roots are —a and — 1(1 ± y 1 — ^|). 

The system (M[) is well-posed if and only if for |£| tending to +oo the real part of the 
eigenvalues associated to A(0 stay non positive. Hence, we must have: 

e, v, a > 0. 

Let us now state a necessary and sufficient condition for the global stability of (M'). 

Proposition 4.8 The linear system (M') is globally stable if and only if the following 
conditions are verified: 



(*) 



u,e,a>0, a(3>0, jS(u + a) + v(3 > 0, 7* + > 0. 
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If all the inequalities are strict, the solutions tend to in the sense of distributions and 
the three eigenvalues Ai(£), A+(£), A_(£) have the following asymptotic behavior when £ 
tends to 0: 



Proof : 

We already know that the system is well-posed if and only if v, a > 0. We want that all 
the eigenvalues have a negative real part for all £. 

We have to distinguish two cases: either all the eigenvalues are real or there are two 
complex conjugated eigenvalues. 

First case: 

The eigenvalues are real. A necessary condition for negativity of the eigenvalues is that 
P(X) > for X > 0. We must have in particular: 



P 5 (0)=ae + ^>0 V^O. 



This imply that aj3 > and as > 0. Hence, given that a > 0, we must have f3 > and 
e > 0. For £ tending to 0, we have: 

|£|2 

Making A tend to infinity, we must have P%{\) > and so j6 + (3 > 0. 
The converse is trivial. 

Second case: 

P^ has two complex roots z± = a ± ib and one real root A, we have: 

Pt(X) = {X - \){X 2 - 2aX + \z±\ 2 ). 

A necessary condition to have the real parts negative is in the same way that P%(X) > 
for all X > 0. 

If 7<5 + (3 > 0, we are in the case where £ tends to (and we see that P^ is increasing). 
We can observe the terms of degree 2 and we get: A + 2a = —a — v then A and a 
are non positive if and only if P^(—a — v) < (for this it suffices to rewrite P^ like 
P ( (X) = (X - X)(X 2 - 2aX + \z±\ 2 )). Calculate: 

2 u(5 + uj5 + ajS 
P({-a — v) = -ve -v a ^ • 

With the hypothesis that we have made, we deduce that P^{—a — v) < for £ tending 
to if and only if v(3 + v^b + aj5 > 0. 
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Behavior of the eigenvalues in low frequencies: 

Let us now study the asymptotic behavior of the eigenvalues when £ tends to and all 
the inequalities in (A) are strict. 

We remark straight away that the condition j6 + (3 > ensures the strict monotonicity 
of the function: A —* Pf(X) for £ small. Then there's only one real eigenvalue Ai(£) and 
two complex eigenvalues A±(£) = a(£) ± 

Let e~ < < e + <0. When £ tends to 0, we have: 

n(A) ~ ier 2 (A( 7 <5 +/?)+<*/?)). 

Then P^(e _ ) < and P^(e + ) > and has a unique real root included between e~ and 
e + . These considerations give the asymptotic value of Ai(£). 
Finally, we have: 

\t(0 + 2a(0 = -a - i/ and - (a(£) 2 + 6(0 2 )A(0 = < + j£ ~ £p 
whence the result. □ 



We summarize this results in the following remark. 

Remark 2 According to the analysis made in proposition we expect the system (M) 
to be locally well-posed close to the equilibrium (p, 0, T) if and only if we have: 

(C) p(p, 9) > 0, \(p,0) + 2fi(p, 9) > 0, K(p) > 0, and X (p,f) > 0. 
By the calculus we have: 

P = dpPo{p,T), 7 = — — =-, 5 = = . 

pd T eo(p,T) p 

We remark that j5 > if dTeo(p,T) > 0. In the case where n verifies rj(p,T) > 0, the 
supplementary condition giving the global stability reduces to: 

(D) d pPo (p,f)>0. 

Now that we know the stability conditions on the coefficients of the system (M ) , we aim 

JV 

at proving estimates m the space E 2 . 

We add a condition in this following proposition compared with the proposition H]8] which 
is: 7$ > 0, but it's not so important because in the system {NHV) we are interested in, 
we have effectively 7<5 = ^77^ > 0. 

Proposition 4.9 : Under the conditions of proposition with strict inequalities and 
with the condition j5 > 0, let (q,d,T) be a solution of the system (M ) on [0, T) with 
initial data (qo,uo,%) such that: 

q 6 P> _1 ' s ,d e P s ~\ T G B 8 - 1 ' 3 ' 2 for some s£l. 



18 



Moreover we suppose that for some 1 < r\ < +00, we have: 
We then have the following estimate for all r± < r < +00; 

+ ||r || 5s _x, S - 2 + ||F|l i?( - s _3 + ^, s - 2+ ^ ) + ll G ll J; r 1(B .-8 + A ) + 11^ 11^-3+^,-4+^ 

Proof: 



We are going to separate the case of the low, medium and high frequencies, particu- 
lary the low and high frequencies which have a different behavior, and depend on the 
indice of Besov space. 



1) Case of low frequencies: 

Let us focus on just the first three equation because the last one is a heat equation that 
we can treat independently. Applying operator A; to the system (Mi), we obtain then 
in setting: 

qi = A t q, d t = A t d, T t = A{T 

the following system: 

d tqi + kdi = F u (4.1) 
d t di - vAdi - eA 3 Ql - pA Ql - 7 A7? = A^divG,, (4.2) 
dfli - aAT t + SAd l = H h (4.3) 

Throughout the proof, we assume that 5^0: if not we have just a heat equation on (|4.3p 
and we can use the proposition 13.51 to have the estimate on T and we have just to deal 
with the first two equations. Denoting by W(t) the semi-group associated to (|4.1l - nO|) 
we have: 

lit) \ ( q \ t ( F(s) \ 

u(t) = W(t) u + / W(t-s)\ G(s) ds . 
9(t) J V #0 / Jo \ H(s) J 

We set: 

fi = PMb + H\\b + ^WiWh - 2K(A<n,di) 

for some K > to be fixed hereafter and (•, •) noting the L 2 inner product. 
To begin with, we consider the case where F = G = H = 0. 

Then we take the inner product of (14, 2p with di, of (|4.ip with f3qi and of (14. 3D with 7^. 
We get: 

IjtWh + m\\b + %Pi\\b) + - e(A 3 q h d t ) + ^WVMb = 0. (4.4) 
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Next, we apply the operator A to (14. 2h and take the inner product with qi, and we take 
the scalar product of (|4.1|) with Adi to control the term 4r(Aqi,di). Summing the two 
resulting equalities, we get: 

^{Mi,di) + [|M||£ a - i/(Ad,,Ag,) -ellA 3 ®!^ - (3\\A<li\\h ~ 7<A^, Am) = 0. (4.5) 
We obtain then in summing (|4.4j) and f|4.5|) : 

5^/? + (HlV^llia - + (KP\\Aqi\\h + Ke\\^1i\\l>) + ^WMb 

zat o uq\ 

+ Ku(Adi,Aqi) + K~j(AT h Aq t ) - e(A 3 q h d t ) = . 

Like indicated, we are going to focus on low frequencies so assume that I < Iq for some 
la to be fixed hereafter. We have then Vc, b, d > : 



|(Ad,,Ag,)| < ^||A^||2 2 + 1|| Ad z ||2 2 



(4.7) 



f72 2 ' r 

|<A 3 «,di>| = KA^M)! < -^"IIAffllli, + gllMlli- 



Moreover we have: ||Vdj|| 2 2 = HMH^a- Finally we obtain: 

+ K\j3 + eC2 2 ^ - u b -- £ ^ - 7^] \\A qi f L2 < 0. 
Then we choose (b, c, d) such that: 

2v e 27 

which is possible if 7 > as v > 0, e > . In the case where 7 < 0, we recall that 7 
and 5 have the same sign, we have then no problem because with our choice the first and 
third following inequalities will be satisfied and if 7 < in the second equation the term 
7^ is positive in taking d > 0. So we assume from now on that 7 > and so with this 
choice, we want that: 

2 

6 f 2 

- +eC2 2l ° -C2 2l ° — > 0, 
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We recall that in your case v > 0, /3 > 0, a > and 7 > 0, 5 > 0. So it suffices to choose 
K and Zo such that: 

2a\ , 2/ / Bu 1 

K < mm 5- , — and 2 < mm I — 



2 (i + C2 2 'o|i)' 5 y V6Ce 2 '6eC 
Finally we conclude in using Proposition 13,11 part (ii) with a c small enough. We get: 

I|/f + c ' 2 2 </f <0 for Z<Z . (4.8) 
2) Case of high frequencies: 

We are going to work with / > l\ where we will determine l\ hereafter. We set then: 

ff = eB\\k qi f L2 + BWdiWl, + HA-^II^ - 2K(Aq l ,d l ), 

and we choose B and K later on. 

Then we take the inner product of (|4.2|) with d\\ 

lj t \\di\\l2+HVd l \\ 2 L2 -e(A 3 qi ,d l )-P{A qi ,d l )- 1 (AT l ,d l ) = 0. (4.9) 

Moreover we have in taking the scalar product of (14. ip with A 2 t#: 

i|||A^|l! 2 + (A 2 ^,A«) = 0. (4.10) 

And in the same way with (|4.3p . we have: 

I^IIA-^Hl, +a||7I||i 2 +5(d l ,A- 1 T l ) = 0. (4.11) 

After we sum (jiH . ((4TTD1I and (|4TTT]l to get: 

I|(£||^||| 2 + e£||A^|li 2 + IIA-^Hl,) + Bi/UVdjUia + ami, 



(4.12) 



S/3(Ag z , d,) - B 1 {AT h di) + 5(di,A~ 1 Ti) = 0. 



Then like previously we can play with (Aqi,di) to obtain a term in ||A^|| 2 2 - We have 
then again the following equation: 

j t (kqiA) + WiWh -v(&d h A qi ) - e\\A 2 qi \\l 2 - f3\\AQi\\h -<y(AT h Aqi) = 0- (4.13) 
We sum all these expressions and get: 

~Ji + [*"l|V4[|i» - A-||M[li.] + mWh + K[fi\\Am\\h + 4^ 2 Qi\\h] , 

1 at - (4.14) 
- B(3{A qi A) ~ B 1 {AT U di) + S(di,A^ 1 Ti) + Ku(Ad h A qi ) + ~fK(AT h A qi ) = 0. 
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The main term in high frequencies will be: ||A^(/i|lx,2- The other terms may be treated 
by mean of Young's inequality: 

KA9l,*)l<^l|A«l| 2 I2 + fl|Adi||| 2 , 

S^IA'Blli. + fllMfe. 
We do as before with the others terms in the second line of (I4.14|) and we obtain: 

\j t ff + (Bu - K)\m\\h] + «\\Ti\\b + K ( JzT + z)W A2( li\\h < 

MiWMh + +K[^\\A 2 gi \\h + y HMIIi, + ^wnb + ^||A 2 *|l! 2 ] 

We obtain then for some a, b, c , d, e to be chosen: 

11 ff + [Bis - (K + B~,\ + + B ^ + ,£-^-)] 

+ [Q -( B7 i- +7 ^ + ^-2_)]| R r ia (4.15) 

+[H+^-^4-74-^i]ii^ii 2 ,^o. 

We claim that a, 6, c , d, e, l±, K may be chosen so that: 

^ _ (# + B7 | + ^2 + B0 d + 5 |_L) > o, (4.16) 

«-(^ + 7^ + *l^)>0, (4.17) 

H+^-^4-74-^^>0. (4.18) 
We want at once that for (|4.16l) and (14.181) : 

V ~l\-P\ >0, (4.19) 
b c 

^-^2"^2 >0 ' (420) 



So we take: 



e = 1 a = 2h6, d = 2h, h = —— — , b = 28h 

c = 2b{y + a)h and h 



2(7(5(i/ + a) + vp) 



With this choice, we get (|4.19p and (|4.2U|) . In what follows it suffices to choice B,K 
small enough and l\ large enough. We have then: 

/, ~ Max(l,2')||^|b + \\&i\\l* +Min(l 5 2')||7I|| L2 
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We have so obtain for I < Iq, I > l\ and for a c small enough: 



I^. / 2 + c ' 2 ^ / 2 <a 

2dt H n ~ 

3) Case of Medium frequencies: 

For Iq < I < li, there is only a finite number of terms to treat. So it suffices to find a C 
such that for all these terms: 



(B) 



IqiWlul 2 ) < C, 



L r T[L2) < C, ||^|| L r (£a) < C for all T G [0, +00] 

and r £ [1, +00] 



with C large enough independent of T. 

And this is true because the system is globally stable: indeed according to proposition 
14.81 we have: 













< p-ci(£)t 


: 


Va, b,c £ L 2 




L 2 




L 2 



withci(^) = min 2 ; 0< ^| <2 ; 1 (Re(Ai(^)), Re(A2(£)), Re(A3(£))) where the Aj (£) correspond 
to the eigenvalues of the system. We have then in using the estimate in low and high 
frequencies in part I4.ll and the continuity of ci(£) the fact that there exits c\ such that: 



So that we have: 



IftWIIra \ \\ 
h Wills I dt I < 



And so we have the result (B). 



ci(0 > ci > 0. 



\\(Qo)i\\l* 
\\( u o)i\\l 2 
\\(%)i\\l 2 



for Iq < I < l\. 



4) Conclusion: 

In using Duhamel formula for W and in taking C large enough we have for all I: 
max(l,2')|| % (t)||L 2 + N/(*)ll^+min(l,2-')||7I(t)|| i2 < Ce"^' ( max(l, 2')||(<7o)* 11^ 

c22! ('-)(ma X (l,2 2i )||F i || i2 + ||Q|| L2 
+ min(l,2- l )||fr,|| La )ds. 



||(d ) i || L 2+min(l,2- i )||(T ) i || L 2)+C f e 



Now we take the L r norm in time and we sum in multiplying by 2 l ^ s for the low 

frequencies and we sum in multiplying by 2^ s+ ~^ for the high frequencies. 
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This yields: 



+ \\do\\ B .-i+j2* 



cT / r-T 



(--1+2) 



KO 



VO 



— ll?0 || J^s — l,s + H^O || l,s-2 

e c ^)(||^(r)|| L2 + HGKr)!^ + WH^W^drY dt 



l>0 



o vo 



e^dlV^^II^ + IIGK^II^ + HA- 1 ^^)!!^)^ dt) . 



Bounding the right hand-side may be done by taking advantage of convolution inequal- 
ities. To complete the proof of proposition 14.9} it suffices to use that u = — A _1 V(i — 
A~ 1 divO and to apply proposition 13.51 □ 

4.2 Global existence for temperature independent coefficients 

This section is devoted to the proof of theorem 12.11 and 12.31 Let us first recall the spaces 
in which we work with for the theorem I2JJ 



E s =[C b {W + ,B s - l ' s )C\L l {m + B s+l > s+2 )} x [C b (R + ,B s - 1 ) N r\L 1 (R + ,B s+1 ) N ] 
x [C b (R+, B 8 - 1 ' 3 - 2 ) n L 1 ^, B s+1 ' s )}. 

In what follows, we assume that N > 3. 
Proof of theorem I2.lt 

We shall use a contracting mapping argument for the function ip defined as follows: 



( F(q,u,T) 

Mq,u,T) = W(t,-)* | n ) + / W(t-s)\ G(q,u,T) | da. 

% J Jo V H(q,u,T) 

In what follows we set: 

p = p(l + q) , 8 = 9 + T ,f= ^(O). 
The non linear terms F,G,H are defined as follows: 
F = — div(gn), 



(4.21) 



G = - u .Vu + V(^\Vp\ 2 ) + 



p P 



An + 



C(p) C(p) 



Vdivu 



+ (V((K(p)-K(p))Ap) + 
Pi(p) PxiP) 



P' Q (p)+TP[{p) P^(p)+fP[(p) 



P 



(4.22) 



Vp 



+ 



v0 , A , (p)VpdivM (du + Vu)p'(p)Vp 
P P 
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where we note: £ = A + fj,, and: 



VP P 



TP 1 (p) TP 1 (p) 



divu-u*.V9+ D ' VU . (4.23) 



1) First step, uniform bounds: 

Let: 

V = \\qo\\ s »-i,N + KII B f-i + \\%\\ s f-i,» ■ 

We are going to show that ip maps the ball B(0,R) into itself if R is small enough. 
According to proposition 14.91 we have: 



\\W(t,-)* n ML k <Cr,. 



E~2~ 



(4.24) 



We have then according (I4.2ip . proposition 14.91 and 14.241 

U(q,u,T)\\ E N <C V +\\F(q,u,r)\\ L1{Sf _ ltf) + \\G(q.u/T)\\ iiin , 



L 1 (_B"2" _ ) 



+ \\H(q,u,T)\\„„ lt 



(4.25) 



Moreover we suppose for the moment that: 

(W) I|9|Il°°(MxI w ) < I/ 2 • 

We will use the different theorems on the paradifferential calculus to obtain estimates on 
\\F(q,u,T)\\ L s *_ lt N , \\G(q,u,T)\\ rim N_^ and {{H^^T^^^n^ 



L 1 (_B"2" _ ) 



1) Let us first estimate ||F(g,n,T)|| ~at_ 1 jv . According to proposition 13.61 we have: 

L yB ~2~ ' 2 1 



and: 



||div(ou)|| , ~n , jv < ||<m|| , jv + ||gu|| , jv ,, 



loull jv < IIqII „ jv llull „ jv 



v , ,. \\U\ 



k/n iv, 1s < o , jv I « . jv,, + 

CTiV JV i -i 7-»— t^— — -!- ] n — 4-1 fp 

Because -B 2 > 2 "t" 1 ^> B 2 and £> 2 > 2 +i ^B2 +1 (trom proposition 13. 2|) . we get: 



||di V M|| Ll(gf _ h N ) < 



_JV JV 11 , ki , , JV , -,, + 
l L°°(B^"'T+ 1 ) 11 H Li(BT +1 ) 



_JV JV li „ JV . 

l L 2 (Bir'ir+ 1 ) li n L 2 (Bir) 



2) We have to estimate \\G(q, u, T)[| jv_ 1 . We see straight away that: 

L (B 2 J 



r Mp) _ P(p) 
P P 



An = K(q)Au 
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for some smooth function K such that K(Q) = 0. Hence by propositions 13.71 13-61 and EO 
yield: 



-Kp) p(p) 



}Au 



< 



Z,i(BT _:L ) 



\\K(q)\\ 



JV . \\U\ 



<l 



<l 



~iV , JV \\U\\ , JV 

'"2") "L^B^"" 1 " 1 ) 



In the same way we have: 

„ r C(p) C(p) 



]Vdivw|| 



< 



~N i JV \\U\\ , , £m, , 

L°°(B"2" '"2") ll L 1 (_B^" +1 ) 



||V(K(p)-^(p))Ag)|| ii(Bf _ 1) < 



l L°°(S^") lnlll L 1 (B^" + ^) 



+ 2n> 



r ^o(p) W) 



]Vp|| l1(b4 _ 1} < 



< 



~JV -, AT , -AT,, JV . 

>L CX '(B^'~ 1 '^T) L 1 (B^ + 1 '^' +2 ) 



After it remains two terms to treat: 

-fp[( P ) fp[(p-y 



p 



Vp|| 



< 



■Flip) FM 



+ 



^l(p) 



)0]Vg|| 



P P 

^(p) ^i(p) 



L 1 ( J B'T- 1 ) 



P 



) V(? iIl 1 (^- 1 )' 



^l(p) _M(p) 
P P 



Vp|| 



L 1 (S"2" _1 ) 



< 



<L a °(B~?- 1 y Hl 'L 1 (B~T +1 ) 



L^B-T- 1 ) 



+ \\Fi( q )TV q \\ Li(Bf _ iy 



According to proposition 13.71 we have: 



IITVoll jv , < ||T|| „ ~jv v , lloll „ ~jv jv , 
11 ^"^(Bir- 1 ) - 11 ll L 2 (BTr'TT- 1 ) ll ^ ll L 2 (Bir-ir+ 1 ) ' 

ll^(g)rv g || zl(B ,_ 1) <a||,|| ioo(s43 nrvg|| L1(B4 _ 1) 



Therefore: 
u[ TP[(p) TP[(p) 



]Vp|| 



JV < I 

L 1 (_B"2" _ ) ~ 1 
+ (1 + 



In the same spirit: 

,.P x {p) Px(p) 



l Loo( ~ f _ lif) | m i Li(sf+1) 



■ ' JV 1 A" J I ^ I I , \ . \ | II I / II - . \ . V i j 

L°°(B^"~ 1 '^") y " ll L 2 (B^"'^"- 1 ) ll ^ ll L 2 (_B^"'^"+ 1 ) 



.)imi 



p 



)V0|| , JV , 



< 



~JV , JV 17" _JV , . JV 
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■ _ N 1 <\\ q \\ N N \\T BF \\ o n 



+ \\q\\ LX{S N- 1 ,N ) \\THF\\ Ll{B » ) , 

where we have: 

T BF = and T HF = ^ A,T . 

KO l>0 
Next we have the following term: 

||n*.Vtf|| , n , < \\u\\ n , lltill , n ,1 . 
And finally we have the terms coming from div(D) which are of the form: 

X'Wpdbn ||L(g)Vpdivu|| * x + ||^Vpdivu|| , Wl 

where we have set: 

L(gi) = A^l + ai)) A'(p)_ 
p(l + xi) p 

Afterwards we can apply proposition 13.71 to get: 



||Vpdiv«|| rl ,_#_i, < ||u 



||L(,)Vpdiv^|| Ll(Bf _ 1} < IlL^ll^^llVpdiHI^^. 
As we assumed that (TC) is satisfied, we have in using proposition 13. 7t 



\\L(q)\\ n <C\\q\\ n . 



So we have: 

A'(p)Vpdivu 



i , n i S, \\u\\ , iv ii ||g|| ~iv , jv (1 + IIqII ~n , at 1 

11 p "L^B^ 1 ) ~ 11 ll L 1 (BT +1 ) ll ^ ll L°°(B^"- 1 'T) v '^'^^(B^"- 1 '^)^ 

In the same way we have in using 13.61 13.71 and 13.21 

(du + Vu)W(p) 

11 p "L^B^"" 1 ) ~~ 11 "L^BT" 1 " 1 )" "L°°(B"2") V ii,n L»(BT)' 



l|V(^|Vp| 2 )|l < \\& - ^)|Vp| 2 )|l , + ||^|Vp| 2 l 



2 ' ri ' "l 1 (b^" -1 ) ~ " y 2 2 " ' ri ' "b 1 (bt2") ' 11 2 1 ' ri "i't^)' 

<ll^)ll L , ( ^ ) lll v H 2 ll, 1( ,f ) + IH v H 2 ll L1( ^ ) - 
<lkll LOO(Bf) l|Vp|| 2 2(sf) + ||Vp|| 2 2(Bf) , 



L°°(BTT) lml L 2 (BTr +1 ) ln " l L 2 (BTT + 1 ) 
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where L{q) = ^ - ^ 



3) Let us finally estimate \\H(q,u,T) 



— N -, N 



\ dMxi p p)Ve) - Ue\\ Li{S ^_ 2) < ( 9 )div(^ l(g )V(?)|| Ll(5 ^_, 4 _ a) 



+ ||div(K 1 (g)V0)|| il( ~ f _^ _ 2) + ||K(g)A0|| Li(i _ f „_ 2) , 



and we have: 



So finally: 

|| div(x(p)w) - *a*|| x , , 2 < ||g|| s 1 * ||T|| i „ 

x (2 + lloll ~at , jv ) . 

Next we have: 



,,,wMp) gPi(p) ^. II < l|/r ,. I, 

+ ||rL 1 (g)divu|| Ll(5 ^_ li# _ 2) + ||^i(9)divu|| Ll( -^_ li4) , 

where we denote: 

,_ P!(p(l+x)) Px(p) 

On one hand, 

II Tdivu II , ~n , at , < ||T|| — n i iv o ||w|| , iv,, , 
||Li(g)divu|| ||Li(g)|| , w , lltill , n , 

whence the desired result: 



g>i(p) ^. || < || || || || 



+ IITII ~jv , at 9 J|«|| jv , 1s (l + llgll ~n , at „ ). 

We proceed in the same way for the others terms which are similar, and we finish with 
the last two following terms: 

II U *- W IIli ( bS-i.S-*) ~ l|T|l Li(Bf ) l|u|l L-(Bf- 1 )' 

U^^U^S*-.*-) £ : Vu|| Ll(5 ^_ li4 _ a) + ||Vu : Vu|| Ll(5 ^_ li4F _ a) . 
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and: 



\\K(q)Vu : Vull , ~jv , jv , < IIqII n \\u\\ „ n . 



so the result: 

„£> : Vit, 



A. < (l+||g|L„ /S #_i,#J||«||.., -V • 



Finally in using (|4.24|) . (|4.25|) and all the previous bound, we get: 

U(q,u,T)\\ E N <C((C+l) V + R) 2 . (4.26) 

Let c be such that || • || n < c implies that: || • ||l°o < 1/3. Then we choose R and rj such 
that: 

fl<ini((3C)-\c,l), and r, < . 
So (W) is verified and we have: 

i/)(B(Q,R)) C 5(0,12) . 

2) Second step: Property of contraction 

We consider (q 1 ,u 1 ,T 1 ), (q 2 , u 2 ,T 2 ) in B(0,R) where we note: 

e l = % + e , T i = y- 1 (e i ) 

and we set: 

(6q = q 2 - q[, 5u = u 2 - u Xl 5T = T 2 - T[) . 
We have according to proposition 14.91 and (|4.2ip : 

){q 2 ,u 2 ,T 2 )- ip {qLjULt r L) (q 1 ,u 1 ,T 1 )\\ e n < 
\\F{q 2 ,u 2 ,T 2 ) - F(qi,ui,Tx)\\ L s N_ lt N 



+ \\G(q 2 ,u 2 ,T 2 ) - G{q 1 ,u 1 ,r l )\\ L 

+ \\H(q 2 ,u 2 ,T 2 ) - H(q 1 ,ui,T 1 )\\ ~N_ t 



(4.27) 



where we have: 

F(q 2 ,u 2 ,T 2 ) - F(q 1 ,m,T 1 ) = -div(q 2 u 2 ) + div(gim) 



G(q 2 ,u 2 ,T 2 ) - G(gi,ui,7i) = 5u*.Vu 2 + uJ.Vtfu + V(\{K' P2 - K' pi )\V p 2 \ 2 ) 



+ V(^<(|V P2 | 2 - |V Pl | 2 ) - p{p)A5u + ^-A5u - C(p) Vdw(5u) 
2 F1 pi 



+ V(if (pi)A(p6q)) + V((^( P2 ) - #( Pl ))Ap 2 ) + _ V5q 

P2 Pi 
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+ - ^)V(ty) + ^V(5g) + TP[(p)V5q + pA5TP[(p 2 ) 

Pi P P 

+ pAe 1 (p' l ( P2 ) - p' l ( P i))Vq 2 + j[Pi(p2) - Pi{ P xWe 2 + j[Pi( P i) - Pi(p)]v6t 

+ pAe 1 P[( Pl )V5q + (^M _ M^l) V/02 divn 2 + (pV^drro, 

P2 Pi Pi 

+ V/9idiv(5u) + - ftl^L) (du 2 + Vu 2 )Vp 2 + ^-^{d(5u) + V5uV Pl 

P2 Pi Pi 

+ pV(5q)(du 2 + Vu 2 )). 
And we have for the part pertaining to H: 

H(q2,U2,T 2 ) - H{q 1 ,u X ,T l ) = (— - — ) div( X (p 2 )V# 2 ) + — div((x(p 2 ) - x(pi))V0 2 ) 

P2 Pi Pi 



+ — div( x (pi)V<5T) + 
Pi 



^i(pi) A(P2) 

Pi P2 



-^-divni H 1 ^ — -divni — 5v* .V62 

A P2 A 



+ ^divcSw - u\V8T + {—- —)D 2 : Vu 2 - —D x : VSu - — (D 2 - D x ) : Vu 2 

P2 A p 2 pi pi pi 

Let us first estimate \\F(q 2 ,U2,T 2 ) — F(qi,ui,Ti)\\ ~n_ 1 n . We have: 
\\F(q 2 ,U2,T 2 ) -i ? (gi,ui,Ti)|| [i( -^_ 1 ^ ) < \\div((q 2 - qi)u 2 )\\ L s % _ x v 

+ ||div(gi(t*a - ui))|| -^.^ , 
< ||^|| i2(i?4) ||n 2 || L2(i?4) + \\Sq\\ Loa{B ^\\u 2 \\ LHB ^ +1) + \\S q \\ L2iB , +1) \\u 2 \\ L2{B ^ 

Next, we have to bound 11^(^,^,7^) — G(gi,iti,Ti)|| ^-:k- We treat only one typical 

L [B~^ ) 

term, the others are of the same form. 

We use essentially the proposition 13.71 to treat the product and the composition, so we 
get : 

uP(Pi) 



A(n 2 - uOll . , < (1 + ki|L , kJ||<HI , v 



1 > • 



Bounding ||/7((7 2 , it 2 , 7^) — H(q%, u%, ~iv_ 1 iv_ 2 is left to the reader. So we get in 

using the proposition 14.91 : 

-9(Qi^X)\\ E ^ <C\\(5q,5u,5T)\\ E N (iKgU;,^')!!^ 

+ ll(Q2>«2>' 7 2)ll E f +2||(gL,ML, / 7L)H i _ ; 4 )• 
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If one chooses R small enough, we end up with in using (14.270 and the previous estimates: 

\\y(q 2 , u' 2 ,T 2 ) ~ *(<?!, «i, Ti)\\ E % < \ IK&Z- Su, 5T)\\ e m . 
We thus have the property of contraction and so by the fixed point theorem, we have 

N 

existence of a solution to (NHV). Indeed we can see easily that E~ is a Banach space. 
3) Uniqueness of the solution: 

The proof is similar to the proof of contraction, hence we will have the same type of 
estimates. So consider two solutions in E~: (qi,u\,T{) and (q 2 ,u 2 ,T 2 ) of the system 
(NHV) with the same initial data. With no loss of generality, one can assume that 
(gi,ui,7|) is the solution found in the previous section. 
We thus have: 

(^0 lkl|li°°([0,T]xRW) < -• 

Let T be the largest time such that q 2 verifies (TC). By continuity, we have < T < T. 
Next we see that: 

Sq = qi — q\, Su = u 2 — u\, ST = T 2 -Ti 

verifies the system: 
' d t 5q + divSu = F(q 2 ,u 2 ,T 2 ) - F(q 1 ,u 1 ,T{), 

d t 8u - ^A5u - ivdiv5u - pKVA5q + (PUp) + fp[(p))V8q + Pl ,?l V5T 
P P (4 1 (T) 

= G(q 2 ,u 2 ,T 2 ) - G(gi,ui,Ti), 
d t Sr - j;A5T + ^^-dW5u = H(q 2 ,u 2 ,T 2 ) - H(qi,u\,Ti). 

We apply the proposition 14.91 on [0,T\] with < T\ < T and we have: 

\\(Sq,5u,ST)\\^ < AiTJWQqMiST)^ 

where we have for T\ enough small A(T\) < ^. 

And we thus have: Sq = 0, Su = 0, ST = on [0, T{\ for T\ small enough and we 
conclude after by connectivity. □ 

We treat now the specific case of N = 2, where we need more regularity for the ini- 
tial data because we cannot use the proposition 13.61 in the case N = 2 with the previous 
initial data. Indeed we cannot treat some non-linear terms such as HTdivuH^!^^.^ or 
1 1 xz* . 1 1 (i?o- — i ) because if we want to use proposition [3]6l we are in the case s\ + s 2 = 0. 
This is the reason why more regularity is required. 
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We recall the space in which we are working: 

E' = [C b {R + ,B°' 1+£ ')nL 1 {R + ,B 2 ' 3+£ ')} x [C b (R +1 B°' £ ') N nL^m+^^'f] 
x [C b (R+, B°- 1+£ ' ) n L 1 (R+, 5 2 ' 1+£ ' )] 

with e > 0, E' being the space in which we have a solution . And E corresponds to the 
space where we show the uniqueness of solution. 

E' = [C b {R + ,B ' 1+£ ')nL 2 {R + ,B 1 ' 2+£ ')} x [C b {R + ,B°' £ ') N r\L 2 (R + ,B 1 ' 1+£ ') N ] 
x [C b (R + ,B°>- 1+£ ')nL 2 (R + ,B 1 > £ ')}. 

Proof of theorem 12.21 

The proof is similar to the previous one except that we have changed the functional 
space, in which the fixed point theorem is applied. So we want verify that the function 
ip is contracting to apply the fixed point. We denote by (qlj ul-, Tl) the solution of the 
linear system (M ) with F = G = H = and with initial data (qo,uo,To) 
Arguing as before, we get: 

\mq,u,T)\\ E , < C V +\\F(q,u,T)\\ So>1+e , 

(4.28) 

+ \\G(q + \\H(q 

if: 

ko|| So , 1+e ' + ||«o|| 5 o.«' + ll T o|l 5 o,-i +£ ' < V- 

Let us estimate \\F(q,u,T)\\ Ll( ~ 01+s > y \\G(q, u, T)\\ Ll ^ and \\H(q,u,T)\\ Ll{ ~ _ 1+s , y 

we just give two examples of estimates in the space E ,the other estimates are left to the 
reader. 

l|div(gn)|| Ll( ~ 01+£ / ) < ||H| L i (B i) + IIHI L i (B 2+ e ')> 

and: 

II^IIl 1 ^ 1 ) < IklU 2 ^ 1 )!!*"!!^^ 1 )) 



\qu\ 



L i (B2+e ') £ IMIi-wlMlii^W) + llal| i2(B2+E ' ) l|u|k2 (B1) . 



We do similarly for \\G(q, u, T)\\ L1 ,~ e /. . The new difficulty appears on the last term 
\\H(q,u,T)\\ Ll ~ _ 1+s r . In fact it's only for this term that that additional regularity is 
needed. Proposition 13.61 enables us to write: 

ll^div«|| i . 1( ~ 0i _ 1+£ / ) < Wr\\ LOO{SOi _ 1+e/) \\u\\ LHBa) . 

II U *- W Hli(B0,-i+ e ') - II T II l i ( bi,i+ £ ')IMI£°°(B°)- 

To conclude we follow the previous proof. Uniqueness in E goes along the lines of the 
proof of uniqueness in dimension N > 3. □ 
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4.3 Existence of a solution in the general case with small initial data 

In this section we are interested by the general case where all the coefficients depend on 
the density and the temperature except k. In this case to control the non-linear terms we 
need that 9 be bounded, that's why we need to take more regular initial data to preserve 
the L°° bound. 

As the initial data are more regular, we need to obtain new estimates in Besov spaces on 
the linear system (M ). 

Proposition 4.10 Under conditions of proposition \4^S\ with strict inequality, let (q,u,T) 
be a solution of the system (M ) on [0, T) with initial conditions (qo,uo,To) such that: 

q 6 5 s - 1 ' s+1 ,u G B S ^ S ,T £ E? s ~ 1,s . 

Moreover we suppose 1 < r± < +oo and: 

n rVi I ris— 3+— ,S— 1+ — ^ _ / f^S — 3+ — ,S — 2+— x rr _ 7>, I ryS — 3+ — ,S— 2+— 

F£Lji(B r i G£Llji(B r i r i), H£Lji(B r i r i). 

W^e £/ien /iaue i/ie following estimate for all r £ [ri,+oo]: 

H^Z- (B— + ^L^B^+i'S+T) + " T IIl- (5— ^f.'+f) ~ IMI^s-M+i + ||t*o||g.-i,. 

+ ||T ||^_ M + \\F\\^„,_ 3+ ± S _ 1+A) + ||G|| i?(gs _ 3+ 4, s - 2+ 4 ) + \\H\\ Lr ^ a _ 3+A ^ +A) . 
Proof: 

The proof is similar to that of proposition 14.91 Low frequencies are treated as in propo- 
sition 14.91 because we don't change the regularity index for the low frequencies. On the 
other hand in the case of high frequencies the regularity index has changed so that we 
have to see what is new. For the medium frequencies we can proceed as in proposition 

EH 



Case of high frequencies: 

We are going to work with I > l\ where we will determine l\ hereafter. We set: 

ft = eB\\k qi \\l 2 + BWdtWl, + - 2K(Aq l ,d l ) 

where B and K will be chosen later on. 

Then we take the scalar product of (|4.3p with 7/, we get: 

— [^[113+ a[|V7l||2 a +<5(Ad lj 7i>=0. (4.29) 
After we sum flOD, (I4TTU1) and dOS]) to get: 

^(WHia + zBW^lWh + WMh) + BvWVd^l, + a\\VH\\h 



(4.30) 



- Bf3{A qi , di) - B 1 (kT u di) + J(Ad,,7l) = 0. 
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We sum (PDI) and (f4TT3j) and we get: 

~ff + [BuWVd^ -jqMllia] +«I|V^||2 2 + [/3K||Ag z ||| 2 + £j A||A 2 % ||| 2 ] 



2dt 



(4.31) 



5/3(A«, d,) - B 7 (A7I, di) + 5(Ad,, + #i/(Ad,, A®) + 7 ^(A7I, Ag,) = 0. 



We interest us after only to the terms of high frequencies, so arguing as in proposition 
1131 we get: 

Let us assume that: 

»" - (* + + ^ + + if) > 0, 

(1) «-( B ^+7^ + ^)>0. 

We recall that v > 0, and a > 0. Next we want to have: 

b 

£ — V— > 0. 

2 

So we take: & = § (we recall that e > 0). So with this choice we get (1) in taking B, K 
small enough and l\ big enough in following the same type of estimate as in the proof of 
the proposition 14.91 We have then for I < Iq, l>l\ and c small enough: 

Ii./2 + c ' 2 2^2 <a 

2dr z - 

and: 

fi ~max(l, 2')||gj|| £ 2 + ||^|| Z 2 + ||^|| i2 
Next we conclude in a similar way as in proposition 14.91 □ 

In the general case the coefficients depend on the temperature and we have to con- 
trol the norm L°° in order to apply the theorems of composition. This motivates us to 
work in the following spaces: 



FT =[C b (R+,B^- 1 ^ +1 ) nL 1 (M +J Bir+ 1 'T+ 3 )] x [C b (R + , Bi^-^Jn 
L 1 {R + ,B^ +1 '^ +2 )] N x [Cfc(l + ,B^^T)nL 1 (f + ,5^ +1 ^ +2 )]. 
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Proof of theorem 12. '61 



The principle of the proof is similar to the previous one and we use the same notation. 
We define the map ip as before with the same F, G and H except that our coefficients 
depends on the density and the temperature. We will verify only that ip maps a ball 
B(0, R) into itself, the end is left to the reader. 

1) First step, uniform Bounds: 

We set: 

ao = \\qo\\ s N +1 + \\u \\ s f + \\%\\ s N_ lt N . 

We denote {<Il^l^l) the solution of (M') with initial data (qo, uo, Tq). We have so in 
accordance with proposition 14. lUl the following estimates: 

\\(qL,u L ,T L )\\ F x <Ca , (4.33) 



(4.34) 



q,u,T)\\ F N < Ca + \\F(q,u,T)\\ L ^~N_ h N + ^ 

+ \\G( q , u ,T)\\ LUS „ , } + im^m^s*-^ ■ 

Moreover we suppose for the moment that: 

(H) Wq\\l°°{RxR n ) - V 2 an d II^"IIz,°°(RxRJV) - V 2 • 

We will now treat each term: \\F(q,u,T)\\ L x , \\G{q,u,T)\\ ^ -n^n and 

\\H( q ,u,T)\\ LU ^_ hh . 

1) We notice that: 

||div(«n)|| , jv , < HqII o jv \\u\\ „ n , 

< IMI „ « \\u\\ n + ||u|| n ,n \\q\\ N . 

2) After we focus on \\G(q, u,T)\\ ~n_ 1 n ■ We have according to proposition 13.71 
H{p,6) ix(p,0) 

An < \\K(q,T)\\ N u , ~n,i n,, 

p p J L1{S N_ h N - H w> ;ii loo(b ^- } ii n Ll(B ^+i,^ +2) 

5, (IMI AT + 11^11 JV )\\u\\ , -JV,, JV , o . 

~ v "^"l°°(B"3-) " "L°°(BT) ; " ll ii(B-2-+ 1 '-2-+ 2 ) 

We proceed in a similar way for the term: 

(C(p) - C(p))Vdivn. 
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Next we have in using propositions 13.61 and 13.21 



\\ P V(k( P ) - K(p))Aq)\\ < 



~ JV 1 JV _ , I W _, — JV I o JV I o W — JV 1 i> i i 

N N , 1 llgll jv , 2s I (1 + ||q|| ,~jv , jv J 
L°°(S"T- 1 '^" + 1 )"^ II L 1 (B"2" +2 ) / v " "L^fBT^^T+Y 



Next, we have to treat the following terms: 

ll[ Z ^ M -^ a '(p)]Vg|| il( ^_ 1 ^ ) < ||L 1 (g)L 2 (T)V g || ii( ^_ 1 ,, ) 

+ \\Lx( q )Vq\\ Ll{ ~» _ lif j + ||L 2 (T)Vg|| il(5f _ x , f }> 

where Li and L 2 are regular function in the sense of proposition 13.71 And we have: 
\\Li(q)Vq\\ ~n , jv < llgll ~jv , jv llgll jv , 



||L 2 (T)Vg|| , _jv , jv < IITII ~jv . jv llgll 

II ^\ J m^l^Y-l.Y) ~ II ll Z,oo ( - B ^--l,^-)ll1ll L l( B ^-+l) 



Finally: 



TpP^p) _ fpl , < \\ q \\ „ N lN+1 \\T\\ -jv lN \\q\\ , jv +1 



+ (IMI ~^ i + l|T|| ~JV , A | 1 1 -A . \ 



and: 

„, P X {p) P X {p) 



)V0|| rl , 5 * < ||Li(g)VT[l ,,-4^ v + ||L 2 (T)W|I 



>*'(T) pV\T)' '""^(BT-Lf) ~ V * "LifBT-^f) ' »~^J » v |lil(B jY_ li jY ) 

+ ||L 1 (g)L 2 (T)V0|| 



||L 2 (T)V6>|| ~jv ,jv < IITII ~jv ,jv iirii 

||Li(g)V#|| ~jv , jv < ||g|| ~jv , jv,, IITII , -jst n , 9 . 

Finally: 



+ (||T|I ^^.^jv +||g|l oo ^jv_ 1 ,jv +1 J||T 
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After we have the following terms: 

||u*.Vu|| , ~N , AT < ||u|| 2 ~ N Nii • 

And we have the terms coming from div(-D). We will treat this one: 
11 A >- (/ ^ )V/;<llv " ' - „ <||L( g! T)Vpdivu|l + ||Vpdivu|| rl , B N lt 



N 



" p L 1 (B"2"~ '"2" ) ~ 11 ' ' L 1 (_B~2~ ' ~2~ ) 11 r ll X 1 (BTr- J -'TT) 

5, (1 + IMI JV , AT ,-, + IITII JV , AT ) llgll ~JV , N ||u|| , — N r i AT . 

Afterwards in the same way we can treat the terms of the type: 

(du + Vu)Vp // x (p,0) (du + Vu)V9 p' 2 (p,0) \' 2 (p, 0)V6 divu 

1 . cinci 1 . 

p p p 

Finally, we have: 

l|v«|Vp| 2 )|| L1( - 4 _ 1 , 4) < \\{k' p - ^)|vh 2 || l1(54 , 4+1) + IIIVp| 2 || l1( ~ 4 



<(1 + 



'Loo (5 4-i,f+i ) ^i ( /il' 2(5 4 + i,f +2) - 



3) Let us finally estimate \\H(q,u,T)\\ ~n_ 1 n : 

n div(x(p, 60V6 1 ) x(p,^) Afl |, ^ nw vr c tayt/jmi 

II -— A0|| il( ~^_ 1 ^ ) < HK^div^^ r)ve)|| Ll(5 ^_ li4) 

+ ||div(ifi(g,r)v0)|| i . 1( ^_ 1 ^ ) + ||iir(g)A^|| ii(gf 

and in using the propositions 13.71 and 13.61 we get: 



+ (lkll z2(5 f ,f +1) + \\T\\ lH ~n^ +1) )\\T\\ l2{b m +1) . 



Next we have: 

„ f TPy(p) fPi(p) 



) divU ll rl , 5 4-l.^ ~ ll L l(<?) divn ILw 5 4-l, 



A r 



where: 



+ ||Mtf)M r ) di Hl il(5 4-i,f } + ||L 2 (T)divu|| , 



||Li(g)divu|| , ~n , n < llgll ~at . at ||u|| n 
||L2(T)div«|| — at i Af < IITII ~ n , at llull , at .I 
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so we get: 

TP 1 (p) TP 1 (p) 
P P 



)divu|| , ~jv , jv < IITII ~jv , n \\u\\ , jv llgll ~jv , jv 

+ ( ~n , jv + T ~jv -, jv I rU , jv i, 



To end with, we have the last two terms: 

\\U .V0 , ~JV , JV < ll'T'll , ~JV JV jo Hull ~JV , JV , 

" "li^^"- 1 '^-) — " "l 1 (B^" +1 '^" +2 )" "L 00 ^^" -1 '"^) 

..D :Vu.. ., , , . „ . 

< (i + lkll ioo(5 4,4 +1) )lkll L2( ^,4 +1) l|n|| i2(i34+1) . 

Finally we have in using (|4.33|) . (|4.34|) and the previous bounds: 

U(q,u',T')\\ EN/2 <C((C+l)r, + R) 2 (4.35) 

Let c such that || • ||^jv/2 < c implies that: || • \\loo < 1/3 then we choose R and cto such 
that: 

R<mf({3C)-\c,l), a <inf^|. 

So (W) is verified and we have then: 

vp(B(0,R)) C B(0,R) . 

Next one can proceed as in the proof of the theorem l2.ll we have to show the contraction 
of the application tp to use the theorem of the fixed point. 

JV 

The uniqueness of the solution in the space F~ follows the same lines as in theorem 12.11 
The details are left to the reader. 

5 Local theory for large data 

In this part we are interested in results of existence in finite time for general initial data 
with density bounded away from zero. We focus on the case where the coefficients depend 
only on the density with linear specific energy, and next we will treat the general case. 
As a first step, we shall study the linear part of the system {NHV) about non constant 
reference density and temperature, that is: 



(iV) 



dtq + divu = F, 

d t u - div(aVn) - V(6divu) - V(cAg) = G, 
d t T-dbf{(NT) = H, 
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5.1 Study of the linearized equation 

We want to prove a priori estimates for system (JV) with the following hypotheses on 
a, b, c, d: 

< ci < a < Mi < oo, < c 2 < a + b < M 2 < oo, < c 3 < c < M 3 < oo, 
< C4 < d < M± < oo. 

We remark that the last equation is just a heat equation with variable coefficients so that 
one can apply the following proposition proved in |13j . 



Proposition 5.11 LetT solution of the heat equation: 

d t T - div(dVT) = H, 

we have so for all index r such that — y — 1 <r< -y — 1 ^ e following estimate for all 
a G [1, +oo]: 

We are now interested by the first two equations of the system (iV) . 



(N' 



dtq + divn = F 

d t u - div(aVn) - V(6divn) - V(cAq) = G 



where we keep the same hypothesis on a, b and c. We have then the following estimate 
of the solution in the spaces of Chemin-Lerner: 

Proposition 5.12 Let 1 < n < r < +oo , < s < 1, (q ,u ) G B% +s x (B%~ 1+S ) N , 
and (F,G) G L?(BT -2+»+2/n) x (L^(fif -3+ s +2/r 1))W 



Suppose that Va , V6 , Vc belong to L\{B 2 ) and £/ia£ <9iC G Lj,{L°°). 

Lei (g,u) G (L r T (B% + s+2 l r ) n Z|(5T + s+1 )) x ((Z^,(sf +s - 1+2 / r )) JV n (L^^f +s ) jv ) 6e 
a solution of the system (N ). 

Then there exists a constant C depending only on r, r%, A, /i, k, ci, C2, Mi and M2 sitc/i 
that: 

"^'^"z^Cs^-^+V^C 1 - ^HVcll^^)) < IKVgd.tio)!!^^ + ||(v^,G0l| z?(jB ^_3 + . +a/nL) 
+ HVg|| I§?(B4 _ 1+S) ||5 t cLi (ioo) + ll(Vg,n)|| £ , {s4+S) (||Va|| z , + ||V6||^ (B#) 
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Proof: 

Like previously we are going to show estimates on qi and u\. So we apply to the system 
the operator A; , and we have then: 

d t qi + divu t = Fi (5.36) 

dm - div(aVu/) - V(6divn ; ) - V(cA^) =Gi + R t (5.37) 

where we denote: 

Ri = div([a, Ai]Vu) - V([6, Ai]div«j) - V([c, A,]Ag). 
Performing integrations by parts and usinf (|5.36H we have: 

-/ Ul V(cAqi)dx = ~ [ c\V qi \ 2 dx- [ (divui (V«.Vc) + ^-d t c 
j r n 2 at j r n Jm. n 1 

+ c.Vqi.VFi) dx. 

Next, we take the inner product of (|5.37p with ui and we use the previous equality, we 
have then: 

J^rOhllia + f c\Vqi\ 2 dx)+ ( (a\X7 Ul \ 2 + b\div Ul \ 2 )dx = [ ((G, + R^.^dx 
2 at - j r n j rN it" 

+ / ((div«i(Vc.V«) + ^^3 t c + cVg/.VF,)) (fx . 



In order to recover some terms in Aqi we take the inner product of the gradient of (|5.36j) 
with ui, the inner product scalar of (|5.37p with Vg/ and we sum, we obtain then: 



d_ 
~dt 



Vqi. Ul dx+ I c(Aqi) 2 dx= [ ((Gi + Ri).Vqi + \divui\ 2 + u l .VF l 

JR N JR N 

— aVui : V 2 qi — bAqid\vui)dx. 



Let a > small enough. We define: 

kf = ||«i||£a + / {Hc\Vqi\ 2 + 2aVqi.u l )dx . (5.39) 

Jr n 

In using the previous inequality and the fact that a\b\ < \a\ + \b\, we have in summing: 

\t^I + \ f WVutf+ablAq^dx^mhi + WRlh*) 
2 at 2 j r n 

x (a||V<#|| L 2 + \\ui\\ L 2) + ||VFi|| L 2(Q||ti/|| L 2 + ||cV^|| L2 ) + - ||ftc||z,«. ||V^||^ 2 (5.40) 

+ ||Vc||x«||Vffl||L3||Vtii||ra • 



40 



For small enough a, we have according (I5.39|) : 

iU 2 <IKI| 2 +/ Kc\Vqi\ 2 dx<hf . (5.41) 
2 J r n 2 

Hence according to (j5.40p and (|5.4ip : 

\±k? + K2 2l k? < h (HCIU, + \\RtWv + \\VF l \\ L2 )\\d t c\\ L o a \\Vq l \\ L , + 2 2l k?\\Vc\\ L i . 
By integrating with respect to the time, we obtain: 

h(t) <e- K22l %(0) + C f e- K22l{t - T) {\\d t c\\ L ™\\V qi (T)\\ L 2 + WVF^r)^ + ||^(r)|| i2 

J o 

+ ||12 { (r)|U 2 +2 , A:,(r)||Vc(r)|U 2 )dr. 

After convolution inequalities imply that: 

^Iknto,^]) < (2-" A:K0) + (2- 2 ^ 1 + 1 /— Vr-OiKv^^ C,)!!^^^ + 2-^||^|l^(^) 



+ 2 r l|V^ll^(z< 2 )l|di c llz4(i/*>) + l|Vc|U (ioc)||A;/|| L r ([0!T] ) 



Moreover we have: 

C^hK \\V qi \\ L 2 + \\ Ul \\ L 2 <Ck h 

Finally multiplying by 2^"2"~ 1+s+ r)' and using (|5.4ip . we end up with: 

(1-CllVcll^oo))^ ||(VF,G)|| 

||(Vg ,«o)|l B 4_ 1+s + ||Vg|| £??(B ^_ 1+S) ||9 t c|| iMLOo) + £ 2<(f 



(5.42) 



2ez 



Finally applying lemma [T] on the appendix to bound the remainder term completes the 
proof 

y- 2 /(f+s-i)pil <C || a |i w [Mi N + c[|6[| ^.JMU, n 



+ C ll c llz|(Bf)^Hz^(i 3 f-) 



□ 



5.2 Local existence Theorem for temperature independent coefficients 

We recall the space we will work with: 

F T =[dr(BT) n 4(5^ +2 )] x [CrCBT- 1 )^ n 4(B^ +1 ) N ] x [C T {B% " 2 ) n Ljr(-B^)] 
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endowed with the following norm: 
|| {q,u,T)\\ FT = 



L 1 T (B I T+ 2 ) + 



JV + U , JV I i + \\u\\~ N , 



+ \\T\ 



L^(B^T- Z ) 



We will now prove the local existence of a solution for general initial data with a linear 
specific intern energy and coefficients independent of the temperature. The functional 
space we shall work with is larger than previously, the reason why is that the low fre- 
quencies don't play an important role as far as one is interested in local results. 

In what follows, N > 3 is assumed. 

Proof of the theorem 12. 4t 

Let: 

q 



n -q° + <f , P n = p(l + g n ), u n = u° + u n , T n = T° + T n and 9 n = + T n 



where (q° ,u° ,T°) stands for the solution of: 

' d t q° - Aq° = 0, 
< d t u° - Au° = 0, 
d t T° - AT = 0, 

supplemented with initial data: 

9°(0) = q , u°(0) = n , T°(0) = %. 
Let (q n ,u n ,T n ) be the solution of the following system: 
' d t q n +i + div(« n+ i) = F n , 

,fC(p n 



(Ni) 



dtUn-i 



div 



( 



dt^+i - div 



x(p n ) 



T 

1 n 



n+1 



-V 



div(ii" +1 ) )-V(K(p")Aq 



n\ a =n+l\ 



where: 



, (^n+i,-u n+ i,T n+ i) t=0 = (0,0,0), 
F n =- div(q n u n ) - Aq° - div(u°), 

K'n 



G, 



(u n )*.Vu n + V(~-^-\Vp n \ 2 ) - V(^^)divu n + V 



,p{p r - 



X(p n ) 



P" 



A (p n )V/9"divu n + (du n + Vu n ) p (p n )Vp n ^ 



l + q r 



1 + . 



P 

Pi{p n ) 



)divit r 
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+ [Po(p") + T"P l( p")]V," _ Au „ / Kfl j,\ 

+ v ( "'^^y } dtT ( j)) + v (g (p-)A g °), 

H n = v( T — ^).v^ x (p n ) ^div U ™ - kt.w 1 + n " n 

1 + q n p n p n 

-A0 o + div(^^-V0°). 

V 1 + q n ' 

1) First Step , Uniform Bound 

Let £ be a small parameter and choose T small enough so that in using the estimate of 
the heat equation stated in proposition 13.51 we have: 



ll^ll 1 N + \\u°\\ , AT,, , N +2 <E, 

||r°|L n 2 + ||u°||~ n , +\\q°\\~ n <A . 



We are going to show by induction that: 

(Pn) \\{q\u n ,f n )\\F T <e. 

As (^Oj^Oj^o) = (0,0,0) the result is true for n = 0. We suppose now (V n ) true and we 
are going to show (V n+ i). 

JV 

To begin with we are going to show that l + q n is positive. Using the fact that B 2 L°° 
and that we take e small enough, we have for t € [0, T]: 

H? n -<folU«((0,T)xR") £ l|divu n ||^ (B ^ ) + ||div(g"- 1 -u n - 1 )|| L ^ B ^ ) + \\divu°\\ ^n^, 



< 2e+\\q n - 1 u n - 1 \\ , iv +1 . 



and: 

ll^- 1 ^" 1 !! <||g n-1 || , ivJK" 1 !! w iv +1 N + lk n_1 || a , n +1 llu™- 1 !! a/ ™. 

Hence: 

Ik" - 9qIU-((o,t)xrJV) < Ci( 2e + (4) + e)e). 

Finally we thus have: 
111 + Qo\\l°°((o,t)xR n ) ~ Ik™ ~ Qo\\l<*>((o,t)xR n ) < 1 + ? n < ||1 + Qo\\l°°((o,t)xR n ) 

+ Ik™ - Qo\\l°°((0,T)xR n )i 

whence if e is small enough: 

-^<l + g "<l+ l|Po||LO ° . 
2p p 
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In order to bound (q n ,u n ,'T n ) in Ft, we shall use proposition 15,121 For that we must 
check that the different hypotheses of this proposition adapted to our system (Ni) are 
satisfied, so we study the following terms: 

a ~TT^' h -TT^' C - K ^' d ~ TTq™' 

In using (V n ) and by continuity of \i and the fact that /x is positive on [p(l + min(go)) — 
a, p(l + max((7o)) + a] , we have: 

< Cl < a n = < M 1 . 

1 + q n ~ 

We proceed similarly for the others terms. 
Next, notice that: 



||Va n |L v < W^l-n^W N < C\\q n u - . . 



JV 



ivrf-ll^,, < cik"ll z « (B » +1) . 

To end on our hypotheses we have to control dtc n in norm || • ||xi,(i<»)' As -B 2 <^-> L 
it actually suffices to bound ||9fC n || n . We have: 

L T (B~2~) 

d t c n = K'(p n )d t q n = K ' {p n ){d\v(q n - l u n - 1 ) - div(u n )) . 
And we have in using the propositions 13.61 and 13.71 
ll^VXdiv^-V" 1 ) -divfu"))^ f) < ||JiC / CP B )div(9 n - 1 u»- 1 )|| il ^ 



+ ||lf / (p B )div(u»)|| tf)> 



< (i + kl^^XIKII^^ + II^^cb^,: 



+ l|{?n ~ 1|l L|(B4+ 1 ) l|un " 1|l ^(B4) ) - 

We now use proposition 15.111 to get the bound on T n , so we obtain in taking r = y — 2: 

l|Tn|l ^( B 4)nL-( B ^^) ~ (H^H^(B^- a ) + l|V( ^ )ll L §?(B ^- 1 ) (5 _ 43) 

x II^HI 1 , tfj- 
11 "ii,(BT)^ 



JV 
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So we need to bound d n in L^(B 



\\Vd n \\ n , < C\\q n \\ 

Now we show by induction (7-n+i). Finally, applying the estimates of propositions 15.121 
and l5.1ll we conclude that: 



i ir+V5 „ +1 , f „ + i,i lFr{1 _ c(i| »ii i|<rf+1) + ri^, +>) + \\cx i(B ^ 

+ IKII lp(Bf > + II (B *,)) < IKVi^.GJII,, + ll^ll,, (B f- S) • 



(5.44) 



Bounding the right-hand side may be done by applying propositions 13.61 and 13.71 For 
instance, we have: 

W F n\\ L ^(B^) < \\div(q n u n ) || £ i, (B w/2) + lldW^^^/a) + || Ag°|| L i^ /2) . 

Since: 

we can conclude that: 

11^11^(^/2) < C(A)+£+^) 2 - 

Next we want to control the different terms of G n . According to propositions 13.61 and 
15.121 we have: 



(u n )* .Vu n \\ , N , < ||li n || N -, ||li n || . 



||V( ^|Vp n | 2 )|' , v , < ||Kr/")|V/)„|-|| x +IIIV/),,, „ 

<(i + lklL, "JII<f lu> 



After we have: 



|| V (^l)div« n || , n , <||divu n || , N ||V(^^-^)|| 

< C\\u n \\ , , #+1 J|cf || / v . 

_ 11 "Z4(B^T + 1 )"^ n L5?(B"2-) 

|| V( CCP^) )di ii <c||n™|| 1 iv +1 ||(Z n || * . 

After we study the term coming from div(D): 
|| A'(p n )Vp n div(u n ) || 



ft-u < ||L(? n )Vp n div(« n )|| r 1/D# _ + ||V//Miv(t/ 
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We proceed similarly for the following term: 

(du n + Vu n )p'{p n )Vp n 
l + q n ' 

Next we study the last terms: 



[P^p n )+T n P[(p n )]Vq n 

l + q n "^(B^T" 1 ) ~ 11 KH ' H "^(B^ 1 ) 11 v * ' H "L^B^"- 1 ) 



5,T\\q"\\* , » +I|T"|I , l|5»|l , „„ a + W 



N 



r A(p n ) 



-]W n || , «• , < C(||g n || iv +1)||T"|| , N , 



||div(P^-V«°)|| . n , <(l+\\q n \\ n )\\u 
We proceed similarly with the other terms: 



V 



-An , v(^^div(u )), V(K(p n )Aq^ 
Let us estimate now ||i? n || , n , . We obtain: 



+ C||V(-^-).V0 n || w # 2 

11 l + q n ' Z/j,(BT -2 ) 



S (i + ll«"ll I , (B « ) )ll«"ll 1 , (B 4 +1) ll^"ll t . „*-., 



We have after these last two terms: 

+ ||ir 1 (T»)div„»|| i , (ij „_ v 

with and l^i regular in sense of the proposition 13.71 and: 

\\K{q n )d\vu n \\ , iv 2 <T||g n || * ||u n || * , , 

|iir 1 (r n )divu n |Uw BJ v /2 - 2 x < |ir n || N . \\u n \\ , iv +1 

II IV V IIL T (B / ) ~ II "l^BTT- 2 )" II £*(B1T+ 1 ) 
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so finally: 



r^-Wii , ,,<(tii 



P" 



and, since iV > 3: 



I _i_ NT"!! "\ IL, n ll 

N T" W-l \ N o I \\U \\ N i 



_i_ ll« n ll NT 71 !! 

11 L^?(B~2~) L^? (B~2" ) "z4(B^" +1 ) 



|« n .vr|| , , iv a < ||n n || , n u \\e n \\ . N . 



-TV. 

.1 . V7„ ,1 I 



|Vu n : Vu r ' 



< u 



T 

mi2 



We obtain in using (I5,44|) and the different previous inequalities: 

\\(q n+1 ,u n+1 ,T n+1 )\\ FT (l - C2^(A + y/i)) < d(e{A + V^? + T(A + Je)). 

In taking T and e small enough we have (V n +i), so we have shown by induction that 
(q n ,u n ,T n ) is bounded in Ft- 



Second Step: Convergence of the sequence 

We will show that (g n ,u n ,T n ) is a Cauchy sequence in the Banach space Ft, hence 

converges to some (q,u,T) £ Ft- 

Let: 

5 q n = q n+l - q n , Su n = u n+1 - u n , 5T n = T n+1 - T n . 
The system verified by (5q n , 5u n , 8T n ) reads: 

( 8 t 5q n + dw5u n = F n - F n ^, 

d t 5u n - div(^^V5u n ^-v(^^dw{5u n )^-V(K{p n )A5q n ) = 

G n — G n -i + G n — G n _i, 

d t 5T n - div( ^^VStA =H n - iJ n _! + H n - <_ l5 

k Sq n (0) = , 5u n (0) = , 5T n {0) = 0, 
where we define: 

G' n = -div((^l - ^W +1 ) - V((K( P ^) - K(pn)Aq^) 



P 



P" 



In the same way we have: 



H 



' n = dw(( xipn+ ]\ -M)V^). 

W 1+ n+ l l + Q n> ) 



+ g n+1 1 + q r 
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Applying propositions 15.11 1 15.121 and using (V n ), we get: 

\\(Sq n ,Su n , ST n )\\p T < C(\\F n - F n -i\\ L ^ B N/^ + \\G n - G n _i + G n - G rt _ 1 || L ^( B jv/ 2 -i) 

+ \\H n - H n -i +H n - H n _ 1 \\ L ^ B N/2-2^ ) ), 

And by the same type of estimates as before, we get: 

\\(Sq n ,8u n ,8T n )\\ FT < CV^(l + A o n(6q n -\5u n -\5T n - 1 )\\ FT . 

So in taking e enough small we have that (q n ,u n ,T n ) is Cauchy sequence, so the limit 
(q, u, T) is in Ft and we verify easily that this is a solution of the system. 

Third step: Uniqueness 

Suppose that (qi,ux,Ti) and (q 2 , u 2 , T 2 ) are solutions with the same initial conditions, 
and (q\,ui,Ti) corresponds to the previous solution. 
Assume moreover that we have: 

||gi(t)lk» <u, yt g [o,T]. 

We set then: 

5q = q 2 — qi, 5u = u 2 - ui, 5T = T 2 - T\. 
The triplet (5q, 5u, 5T) satisfies the following system: 

' d t 5q + d\w5u = F 2 -Fx, 
d t 5u - dw(^^-V5u) - V(^^-dw(5u)) - V(K(p 2 )ASq) =G 2 -G l + G > 
d t 5T - div(^^-V5T) = H 2 -H 1 + H, 

1 + 92 



with: 



k 6q(0) = , Su(0) = , ST(0) = 

G'=- divf {f^L - ^l)Vu 2 ) -V((K( P2 ) - K( Pl ))Aq 2 ) 

V P2 Pi J 

C(P2) C(pi) 



v( (! 



P2 Pi 



-)div(m) , 



H' =div((^-*^)W 2 
Vl + q 2 l + qi J 

Let T the largest time such that: ||?2||ioo((o,f )xR N — a - As 12 G C([0, T]; B N / 2 ), we have 
by continuity < T < T. 

We are going to work on the interval [0, T±] with < T\ < T and we use the proposition 
15.121 so we obtain in using the same type of estimates than in the part on the contraction: 

\\(5q,5u,5T)\\N < Z(T)\\(5q,5u,5T)\\M 

p 2 p 2 
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with Z(T) — >t^o 0. 

We have then for T\ small enough: (Sq, Su, 5T) = (0, 0, 0) on [0, Ti] and by connectivity 
we finally conclude that: 

Qi = 92, u\ = u 2 , T\ = T 2 on [0,T]. 

□ 

Proof of the theorem 12.51 



In the special case N = 2, we need to take more regular initial data for the same reasons 
as in theorem !2.2l Indeed some terms like ^(6)divu or u*.V6 can't be controlled without 
more regularity. 

The proof is similar to the previous proof of theorem 12.41 except that we have changed 
the functional space i*r(2), in which the fixed point theorem is going to be applied. As 
we explain above we can use the paraproduct because we have more regularity, so we 
just see the term u*.V9. The other terms and the details are left to the reader. 
We then have: 

U U *- V *lli4(B-i.-i+«') ~ ll T ll^ ( 5o,i +e ')H w lli??(B )- 

□ 



5.3 Local existence theorem in the general case 

Now we suppose that all the coefficients depend on the temperature and on the density, 
and that conditions (C) and (D) are satisfied with strict inequalities. 
One of the problem in the general case is the control of the L°° norm of the temperature 
6 in order to have control on the non linear terms where the physical coefficients appear. 
Indeed in the theorem of composition we need to control the norm L°°. 

N 

So we must impose that 9q is in 2 to hope a L°° control. And in consequence the 
others initial data have to be also more regular. 



Proof of theorem 12. 61 



We proceed exactly like in theorem 12.41 except that we ask more regularity for the initial 
data. We define then: 

q n = q° + q n , u n = u° + u n , T n = T° + T n and 9 n = T n + T 

where (q° ,u° ,T°) stands for the solution of: 

d t q° - Aq° = 0, 
d t u° - Au° = 0, 
d t T° - AT = 0, 

supplemented with initial data: 

q°(0) = q , «°(0) = u , 0°(O) = % 
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Let (q n ,u n ,8 n ) be the solution of the following system: 
8 t q n+1 + div(u n+1 ) = F n , 



d t u n+1 - div O^p- Vu n+1 ^j - V ( C( y ^ div(n" +1 )) - V(K(p n )A<t +1 ) = G ri 
d t T^-dU^^T-+A=H n 



^(q n+ \u n+ \f n+1 ) t=0 = (0,0,0). 
where: 

F n =- dw(q n u n ) - Aq° - div(u°) 



Gn — 



+ V(^f |V„f ) - V(4^)div„" + V(«9Adi™-' 



A , 1 (p n ,g n )Vp n divn n (dn TO + W 1 ) ^Qo", fl ra )Vp 



1 + 9 



A' 2 (p n ,e n )V0 n divu n + (du n + Vu n )n' 2 (p n ,9 n )Ve n 



l + q n 



+ 



[P^p n ) + T n P[(p n )]Vq n 



l + q n 



+ 



l + q n 

Pl(p n ) 
p n ifj'(T n ) 



vr 



_ AM o + div(^V„°) + v(*^ldiv, u o,) +VW W 

i '-pn p ( n n\ T~\n . V7 ? .fi 

^ =v(— -).vr x ( P ",n ^dW 1 - («»)*. V0 n + — 



-Ag + div( x(p "' gn) Vg°). 
1) First Step , Uniform Bound 

Let e be a small positive parameter and choose T small enough so that in using the 
estimate of the proposition 13,51 we have: 



\\r° 



+ \\q c 



< Ac 



3 <£, 



ir°|L jv at 

i ii roo/n ~T" " 1 1 t oo / d T" " - 1 lire 



After we are going to show by induction that: 
(Vn) \\(Qn,u n ,f n )\\ FT < e. 
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As (qo,uo,To) = (0,0,0) the result is true for (Vq). We suppose now (V n ) true and we 
are going to show (V n+ i). 

To begin with we are going to show that 1 + q n is positive. In using the fact that 

N 

B~ L°° and that we can take e enough small, we have for t G [0, Tj: 
||« n -9b[|L«(( ,T)xR^)< C i(H div (^)ll^ (B 4 ) + l|div( ( ?"- 1 ^- 1 )||^ (B4) 



+ l|dSv(«°)|| 



< Cl (2Vfe + \\ q ^-\ lT{B§+1) ), 



and by induction hypothesis (V n -i): 



^n—X„,n—ln ^. frp\\ _n— 1 11 IL.n— 111 

1 M Li,(S^- +1 ) - M L5?(BT) M II L2(s^-+ 1 ) 



i /r^ll 1 II Wn.n— 1|| 

thus: 

Ik™ - 9o||l~((o,t)xR") < CiVT(2e + (4 + e)e). 

Finally we have: 

111 + 9o||l°°((o,:t)xr jv ) ~~ \\l n ~ Qo\\l^((o,t)xr n ) < 1 + <? n < P + Qo\\l°°((o,t)xm. n ) 

+ \\q n - <?o||l°°((0,T)xIR jv )> 

C / 1 , „n / 1 , llA)l|i° 



< 1 + g n < 1 + 



2p" - p 

So we have shown that: 

(*) \\q n \\L°° < 2A 

and that p n is bounded away from 0. 

To verify the uniform bound we use the propositions 15.111 and 15.121 For that we have 
to verify the different hypotheses of these propositions, so that we study the following 
terms: 

I + q n ' l + q n ' yHnJ l + q n 

In using (V n ) and by continuity of fi and the fact that fi is positive on [p(l + min(q )) — 
a,p(l + max(<7 )) + a] x [8(1 + min(T )) — a, p(l + max(T )) + a], we have: 

0<Cl < o » = 4^<M 1 . 

1 + q n ~ 
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We proceed similarly to verify the bounds of the other terms. 

After we use the proposition 13.61 and the fact that q n is bounded. We get: 



Next we want to estimate dtc n in Lj,(B~). For that, we use the fact that: 
d t c n = K\p n )d t q n = ^'(p")(div(g n - 1 u n - 1 ) - div(u n )) 

And we have: 

||KV)(div(<f- V- 1 ) - div(«*))|| ^ } < llAVjdiv^-V- 1 )!! £ ) 



+ ||if(^)drv(«»)||^ )> 



<(l + ||? n || at )(||u n || , n +1 +|k n - 1 n n " 1 || , N+1 ) 



<(i + |k n || jv XIKH , N+1 +\\q n - 1 \\ n IK" 1 )) , N+1 

~ X llH "Lf?(BTF) nl1 ll L^(B^ + 1 ) " H ll L^(B^) U ll L^(B^ +1 ) 



+ii'"-\ 3 W""~ 1 W> ) - 

Now we want to show (V n +i) by induction and in this goal we will apply the estimates 
of proposition 15.111 and proposition 15.121 This is possible as we have verified above the 
validity of the hypotheses. We obtain: 



+ IKIU (rf) + IIVIL < l|(VF„,G„)|| . + \\H n \\ i} . 



(5.45) 



We want to control now the part on the right-hand side of (|5.45|) . for this we do like 
previously in using proposition 13.61 We have: 

with: 

\\ un Q n \\ i, N +2s < \\Q n \\L° c (L°°)\\u n \\ JV+o, + Ik"!! 0/ , 2 ^ 11^^ II iV N 



< \\q n \\ ~n N+1 \\u n \\ . N+2 +VT\\q n \\~ 9 ~n +1 n, 2 \\u n \\ n . 



One ends up with: 



\\F n \\ LlT{B N +1) <C(A + e + V~e) 2 . 
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Next we want to control the different terms of G n . We have: 

l|v(^|v,„p)|| L , (flf) < ^lk»|| t?(B » +1) )lk ll ^ (fl „« ) . 



After we have: 



We treat similarly the term: 



Next we study the term: 

" l + q n "l^B 4 ?") ~ V " ' 114 11 ^(B^)' V " 11 ** ll 4(BT+ 1 ) ll ' J ll L5f(S'T +1 ) ' 

We proceed similarly for the following term: 

(du n + Vu n )fi'(p n ,e n )v P n 

l + q n ' 

Next we have: 

[Po(p")+f"PV)]WV < 



+ T\\T n \\ , N \\q n \\ , (l+||(? n || , ivj , 



Pl ^ n) -]Vr|| iv < (\fT{\\q n \\ , N +\\T n \\ . n. 



p n ^'(T n ) "4-( BT ) ~ V L^?(B~2~) 11 "L?(B^y 



We proceed similarly with the other terms: 



-An , v(^^div(n )), V( K (p n )Aq°). 
After we want to estimate the term \\H n \\ , n . So we have: 

l|V(:r^— ).w n x(p n ,0 n ))ll w n <(i+lk n ll , ^ + 11^1 , *J 



x VT||(7 n || , iv +1 J|T n || 2/ n 
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We have after these last terms: 
,,T n P x (p 



divu n \\ , n < Vf(\\q n \\ n +\\T n \\ n (l+||g n || N )) 



x ||if n || , , N , , . 



and: 



IK.wi , n < >/f||u n || iv ||r n || , iv +1 

\\Vu n : Vu n || , n < \\u n \\ 2 N 

we obtain in using (|5.45[> . the hypothesis of recurrence to the state n and the previous 
inequalities: 

\\{q n+1 ,u n+1 ,T n+1 )\\ FT (l - C2^(A + V^)< CMA + VI) 2 + T(A + y/e)). 

In taking T and e small enough we have (V n +i), so (q n ,u n ,T n ) is bounded in Ft- To 
conclude we proceed like in the proof of theorem 12,41 and we show in the same way that 
(if 1 ,u n ,T n ) is a Cauchy sequence in Ft, hence converges to some (q,u,T) in Ft- We 
verify after that (p, u, 0) is a solution of the system. 



Uniqueness: 

We compare the difference between two solutions with the same initial data and we use 
essentially the same type of estimates than in the part on contraction. The details are 
left to the reader. 



6 Appendix 

This part consists in one commutator lemma which enables us to conclude in proposition 
15.121 Moreover we give the proof of proposition 13.71 on the composition of function in 
hybrid spaces adapted from Bahouri-Chemin in [2J. 

Lemma 1 Let < s < 1. Suppose that A £ L T (B^ +1 ) and B G L T (B^ _1+s ). Then 
we have the following result : 

||dJAA z ]5||rwr2^ < Cq2- , ($- 1+ ">p|L, n +1 A\B\U n 1+s 

II KL ! tj \\L T (L ) — I II l\ L 2 II ll£,J B^r-H» 



with YjI& c i = L 
Proof: 



We have the following decomposition: 

uv = T u v + T u 
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where: T u v = Si-iuAiv and: T' v u = YLi&l Si+2vAiu. 
We then have: 

dk[A, A{\B = d k T AiB A - dkA^A + [T A , Ai]d k B + Tq^B - A t T dkA B. (6.46) 
From now on, we will denote by (cj)jez a sequence such that: 

£q<l. 

Now we are going to treat each term of (I6.46H . According to the properties of quasi- 
orthogonality and the definition of T we have: 



d k T' AiB A= 9 fc (5 w+2 A z BA m A). 



m>l-2 

Next, in using Bernstein inequalities, we have: 
\\dkT AlB A\\ L ^ L 2) < ^2 2m \\Ai B \\Ll(L°°)\\AmA\\ L ^( L 2) 



m>l-2 

<2'f||A ZJ B||^ (L2) 2- m f(2 m ^\\A m A\\ Ll{L2) ) 

m>l-2 



< 2 -. W2 -i +s)(2 Kf-i +s )|| AzjB | L|(i2)) £ (2-(f+ 1 )||A m ^|| L , (L2) ) 

m>Z-2 

<c,2-^— »ll«ll E , (B »-, + ,MII^ (Bf+1) . 

Next, we will use the classic estimates on the paraproduct to bound the second term of 
the right-hand side of (|6.46p . We obtain then: 

After in using the spectral localization we have: 

\\d k AiT B A\\ L i T ^ L 2) < 2 1 \\AiT b A\\ l i t(L 2) 

< Q 2-^-^)||B|| 4(sf „ 1+3) P|| z , (Bf+1 . 

According to the properties of orthogonality of Littlewood-Payley decomposition we have: 

[T A , A[]d k B = Y, [Sm-iA, Ai]A m d k B . 

\m-l\<A 



In applying Taylor formula, we obtain for x E 



pN. 



[Sm^AAilAmdkBix) = 2~ l f f h{y)(y.S m ^VA(x-2- l Ty))A m d k B(x-2- l y)dTdy . 

JR N JO 
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By an inequality of convolution we have: 

||[S m -iA A z ]A m %B|| i2 < 2-*||V,4|M|A»A£|| L 2 



So we get: 
Finally we have: 



L 2 T (B^~ 1+S ) 



N _ 



Te fcA AiB= S^dkAA^B. 

|«-m|<4 

Hence: 

\\ T d k A^iB\\ L ^ L 2) < \\dkA\\ L ^ Loo) \\AiB\\ L ^ L 2) . 
And the classic estimates on the paraproduct give: 

The proof is complete. □ 
Proof of proposition 13. 71 

To show (i) we use "first linearisation" method introduced by Y.Meyer in [23], which 
amounts to write that: 

F(ui,uz, ■■■ ,Ud) = y^(F(S p+ iUi, • • • , S p+1 u d ) - F(S p ut, • • • , S p u d )). 

According to Taylor formula, we have: 

F(S p+ iUi, ■■■ , S p+1 u d ) - F(S p u 1 , • • • , S p u d ) = m\v{ -\ h rripU p d 

with u P = ApUi and 

m p = diF(S p ui + su\, • • • , SpUi + su P , • • • , S p u d + su p d )ds. 
Jo 



Observe that: 
We have: 



\m l p \\ L °° < ||VF|| Lo 



A p F(u 1 ,u 2 , ■■■ ,u d ) = Ap + Ap 
where we have decomposed the sum into two parts: 

a« = Wa p ( U X) + • • • + ^p(«X)) > 



<?>P 
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Now we bound ||A P \\l p t (lp) m this way: 

ll A p ) ||i4(Li') ^ Yl (W U l\\L p T (LP)\\ ml q \\L^(L^) + ■■■ + ||^||l?,(Lf)II^I|l-(L-)) 

q>p 

< £ l|VF|| L .2-fp J (||ui|| £5 , (B » ) + • • • + ||u d || £ , (B|) ) 
q>p 

with (c g ) G i 1 (Z). 
Therefore, since s > 0: 

E 2P 1A«|| L , < cuvFiu-diu!^^ + • • • + Utidllz^)). 

(2) (2 

To bound ||A^ ; || l p 

(lp) we use the fact that the support of the Fourier transform of A p 
is included in the shell 2 P C, so that according to Bernstein inequality: 

lK 2) || L p (iP) < £ (||A p (uX)ll^(L,) + --- + l|A p (uX)ll^(L,)), 
g<p-l 

< C2 -p(W+D £ (HflW+Hufm;)!!^ {LP) + ■■■ + ||5 W+1 («X)lk(i»)- 
<?<p-i 

Moreover we have according to Faa-di-Bruno formula: 

/•l m 
d fc m* = / Al.. lm F m+ \S q {u) + su«) J] d ln (S q (u) + su«))ds. 

h+-+l m =k,l m ^0 n=l 

Hence we get for all fceN: 

\\9 k m q \\ L ^( L oo^ < C Uit k2 qk 
with: C Uz:k = C(l + |K||lo?(l°°))- 
We have then: 

H 2) H m <C2-*M+V £ ^ 2 «C-^i-W+i) C7ttij ___^ ( || Ul ||_^ (jH5) _, ,_ H^H^^j). 

<?<p-i 

Hence the result: 

^HIAfH^ < ^,.,^(11^11^^ + • • • + IKH Z , (B|) ). 

So the first part of the proof is complete. 
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For proving (ii) we proceed in the same way as before. We get: 

F(u) = y^m q u q . 

And we have for p > 0: 

A p F(u)=Al + Al 

so: 

< E II^IU-2-^c ? ||n|| z , (gn , ii2) . 
Hence in using convolution inequality: 



T 



j2 2PS2 \\ A l\y T m<c\\F'\\ La o\\ u \\ 

LP,(B S 1' S 2)- 

p>0 



After we get for all s > 0: 



\\KH (V) < C2-f(W +1 ) £ \\dW + Hm q u q )\\ LP , L2) , 



q<p-l 



•]+!-»(?)) 



T 

g<p-l 



with s(g) = si or S2- 
So we obtain: 



^2-||A^| L , (i2) <^2-(H^-)^ Cg 2^W^-)||n|| z , ( ^, fl2) 

p>0 p>0 g<0 

, 2 - P ([s]+i-s 2 ) sr^ 2g ([ s ] + i_ S2 )|, _ 
^ 2.^ 4 ii u Hl5,(B'i.»2)- 

p>0 0<<j<p-l 



(6.47) 



We have to choose s, so for the first term of (|6.47p we just need that: [s] + 1 — S2 > and 
[s] + 1 — si > and for the second term of (j6.47H we just have a inequality of convolution. 
So we can take s = 1 + max(si, S2). 

We do the same for p < and we have: 

^2^||A^|| L , (L2) <^2^^||F'||^2-^ C? |K|| £p(gsi , S2) 

p<0 P<0 9>P 

+ ^2^ £ ||F'||^2-^ C(? ||n|| Zp(5si , S2) . 
p<0 p<g<0 

We conclude by a inequality of convolution. 
And for the term we get: 

I> s lAj| b(i2) <^ 2 ~p( W+ i- Si ) £ c^(W^-)||.|| £p( ^, S2) . 

p<0 p<0 9<P-1 
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For proving (iii) and (iv), one just has to use the following identity: 

G(v) - G(u) = (v-u) [ H{u + t(v - u))dr + G' {0)(v - u) 

Jo 

where H(w) = G (w) — G (0), and we conclude by using (i), (ii) and proposition 13.61 □ 



7 Annex: Notations of differential calculus 

If f : M. n -> R , we denote: 

df = y~] -K—dxi = difdxi 

1=1 

with the summation convention on repeated indices and the simplified notation: 

v-SL. 

OXi 

The vector field associated to the differential df is noted V/, 

dxi dxi 

Let / : W 1 — * M. n . Let denote fa the i th component of /, and: 

(df)i,j = djfi. 

By analogy with the case of the scalar, we denote: 

V/ = (df)* , so (V/) M - = difj. 

The curl of f is given by: 

(curl/); j = difj - djfi. 
The divergence of the vector field f is given by: 

div/ = trdf = difi. 
If A : R n -> R nxn , with coefficients a^we set: 

d 

(div A)j = didij , divvl = div(^4— — ) dxj. 

OXj 

In particular, for / scalar, we have: 

div(/J) = df. 

And finally we set: 

A : B = dijbij . 
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